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AFIT/GE/ENG/OlM-19 


ABSTRACT 

A  novel  phugoid  damping  control  design  methodology  is  developed,  based  on  the 
use  of  wind  axes  and  a  point-mass  aircraft  model.  The  state  variables  are  air  speed,  flight 
path  angle,  and  heading  angle,  the  control  variables  are  thrust  setting,  angle  of  attack, 
bank  angle,  and  sideslip  angle,  and  the  command  signals  are  airspeed,  flight  path  angle, 
and  heading  angle  or  heading  rate.  All  the  variables  and  parameters  are 
nondimensionalized.  A  multivariable  set  point  controller  is  developed  which  consists  of; 
(i)  a  trim  calculation-based  nonlinear  feed-forward  control  computer;  thus,  given  a 
commanded  new  trim  state  (air  speed,  flight  path  angle,  and  yaw  rate),  the  required  trim 
thrust  setting  and  trim  angle  of  attack,  bank  angle,  and  sideslip  angle  inputs  are 
determined,  and,  (ii)  a  small  signal  linear  feedback  regulator;  the  equations  of  motion 
linearized  about  the  trim  condition  of  wings  level,  and  constant  altitude  flight,  which 
simplifies  the  dynamics  to  allow  separation  between  the  lateral  and  longitudinal  control 
channels,  are  used,  and  a  small-signal  linear  multivariable  regulator  is  designed.  The 
linear  compensator  also  entails  integral  action.  Thus,  the  controller  consists  of  a  strongly 
nonlinear  feed-forward  module  and  a  linear  small  signal  compensator.  The  novel 
proposed  multivariable  nonlinear  set  point  controller  encompasses  full  three-axes 
autopilot  functions.  Moreover,  this  controller  is  used  as  a  tracking  controller,  a.k.a.  a 
“phugoid  damping”  controller,  provided  that  the  bandwidth  of  the  command  signal  is 
substantially  less  than  the  bandwidth  of  the  closed  loop  flight  control  system.  The 
phugoid  damping  controller’s  performance  is  examined  in  extensive  simulations  and  its 
wide  operational  envelope  is  demonstrated. 
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NOMENCLATURE 


aj  Tail  Lift  Curve  Slope 

aw  Wing  Lift  Curve  Slope 

AR  Aspect  Ratio 

b  Wing  Span 

Cd  Drag  Coefficient 

Cl  Lift  Coefficient 

Cdo  Zero-Lift  Drag  Coefficient 

D  Drag 

e  Error  Signal 

fs  Sampling  Frequency 

F  Side  Force 

g  Acceleration  of  Gravity 

K  Parabolic  Drag  Polar  Constant 

Ki  Controller  Gain  Matrix  -  Integral  Part 

Kp  Controller  Gain  Matrix  -  Proportional  Part 

K. .  Controller  Gain 

L  Lift 

m  Mass 

q  Dynamic  Pressure 

St  Vertical  Tail  Surface  Area 

Sw  Wing  Surface  Area 

t  Time 

ts  Settling  Time 

T  Thrust 

V  Velocity 

z  Integrator  Charge 

a  Angle  of  Attack 

P  Sideslip  Angle 

Y  Flight  Path  Angle 

p  Air  Density 

X  Directional  Control  Channel’s  Time  Constant  (specified) 

p  Thrust  (Normalized) 

(])  Bank  Angle 

(t)p  Disturbance  Bank  Angle 

v|/  Heading  Angle 

m  Yaw  Rate  at  Trim 

“n  Phugoid  Natural  Frequency 

Lateral  Control  Channel’s  Natural  Frequency  (Specified) 

®np  Longitudinal  Control  Channel’s  Natural  Frequency  (Specified) 
^  Phugoid  Damping  Ratio 
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Lateral  Control  Channel’s  Damping  Ratio  (Specified) 
Longitudinal  Control  Channel’s  Damping  Ratio  (Specified) 


PHUGOID  DAMPING  CONTROL 


1.  INTRODUCTION 

Phugoid  damping  control  entails  the  design  of  compensators  for  the  control  of  the 
aircraft’s  slow  states.  Thus,  in  the  pitch  channel,  one  controls  the  airspeed  and  flight  path 
angle.  In  the  directional  channel,  the  heading  angle  is  controlled.  When  one’s  attention 
is  confined  to  the  regulation  function,  one  then  refers  to  autopilots,  viz.,  altitude-hold, 
Mach-hold,  and  heading-hold  autopilots.  Traditionally,  automatic  pilots  have  been  used 
extensively  in  missiles  and  in  aircraft  to  decrease  pilot  workload  and  improve  flight 
safety.  It  is  envisioned  that,  in  the  future,  complex  outer-loop  flight  control  systems  will 
be  used  increasingly  in  the  control  of  emerging  Uninhabited  Combat  Air  Vehicles 
(UCAVs).  Thus,  the  challenge  is  to  increase  the  operational  envelope  of  autopilots,  and 
also  increase  the  bandwidth  of  the  command  signals  that  outer-loop  tracking  controllers, 
a.k.a.  “phugoid  damping”  controllers,  can  successfully  handle,  thus  endowing  UCAVs 
with  high  performance  autonomous  flight  control  systems. 

Model  based  control  design  is  standard  practice  in  aerospace.  Conventional 
autopilot  design,  as  espoused  e.g.,  in  [1],  is  almost  exclusively  based  on  a  linear  plant 
model  derived  from  a  linearization  of  the  6  DOF  (Degrees  Of  Freedom)  aircraft  equations 
of  motion.  One  acknowledges  the  inherently  slow  time  scale  of  the  autopilot  controlled 
variables  by  using  a  nested  control  loop  structure.  The  inner  loops  consist  of  the  SAS 
(Stability  Augmentation  System)  and  CAS  (Control  Augmentation  System),  whereas  the 
autopilot  is  relegated  to  the  outer  loop.  Typically,  control  loops  for  aircraft  control 
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systems  are  designed  from  the  “inside  out”.  That  is,  the  first  step  is  the  design  of  an 
inner-loop  stability  augmentation  system  (SAS),  including  actuator  dynamics.  From 
there,  the  design  moves  progressively  toward  more  complex  outer-loops:  control 
augmentation  system  (CAS),  and  altitude/Mach/heading  -hold  autopilots. 

The  standard  autopilot  is  a  linear  compensator,  augmented  with  gain  scheduling 
[3].  Conventional  autopilots  have  performed  well  over  many  decades.  At  the  same  time, 
the  operational  envelope  of  conventional  autopilots  is  somewhat  limited,  viz.,  the  linear 
control  design  envisages  small  perturbations  in  the  states,  low  amplitude  set  point 
changes,  and  slowly  varying  set  point  settings.  Thus,  in  World  War  H,  the  autopilot 
controlled  VI  cruise  missile  was  downed  by  RAF  pilots  by  exploiting  the  VI  autopilot’s 
limited  operational  envelope.  They  tipped  the  VI  over  with  their  wingtips  [2,10].  Also, 
autopilot  “upsets”  have  been  reported  in  recent  times  [4,6]. 

In  the  present  paper,  the  model-based  design  of  the  autopilot  /  outer  loop 
controller  hinges  on  a  point-mass  model  of  the  aircraft  dynamics.  The  low  order 
dynamics  exclusively  encompass  the  states  relevant  to  autopilot  design.  The  “fast”  inner 
loop  states  are  the  control  variables. 

The  novel  approach  to  outer-loop  controller  design  pursued  in  this  thesis  is  as 
follows:  Directly  design  the  outer  loop  controller  employing  a  low-order,  slow- 
dynamics,  albeit  nonlinear  model.  It  is  hypothesized  that  such  an  approach  would 
simplify  the  outer-loop  controller  design  and  help  better  capture  the  nonlinear 
characteristics  of  the  air  vehicle.  The  latter  is,  in  part,  conducive  to  a  “full  envelope” 
controller,  thus  obviating  the  need  for  gain  scheduling.  Although  this  approach  does  not 
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include  the  “fast”  dynamics  in  the  outer-loop  controller’s  design,  ensuring  that  sufficient 
“phase  margin”  exists  in  the  outer-loop  compensates  for  this  deficiency. 

There  is  reason  to  believe  that  ab  initio  using  a  low  order  plant  model  is  the  right 
approach  to  outer  loop  controller  /  autopilot  design.  Moreover,  using  the  right  plant 
model,  a  nonlinear,  trim  solution-based  controller  naturally  suggests  itself.  Hence,  a  high 
quality  outer  loop  controller  /  autopilot  with  an  expanded  operational  envelope  is 
realized. 

Specifically,  in  this  paper,  a  novel  phugoid  damping  control  design  method  is 
developed,  based  on  the  use  of  wind  axes  and  a  point-mass  aircraft  model.  The  state 
variables  are  air  speed,  flight  path  angle,  and  heading  angle,  the  control  variables  are 
thrust  setting,  angle  of  attack,  bank  angle,  and  sideslip  angle,  and  the  command 
(reference)  signals  are  airspeed,  flight  path  angle,  and  heading  angle  or  heading  rate.  All 
variables  and  parameters  are  nondimensionalized.  A  multivariable  set  point  controller  is 
developed  which  consists  of:  (i)  A  trim  calculation-based  nonlinear  feed-forward  control 
computer;  thus,  given  a  commanded  new  trim  state  (air  speed,  flight  path  angle,  and  yaw 
rate),  the  required  trim  thrust  setting,  angle  of  attack,  roll  angle,  and  sideslip  angle  inputs 
are  determined,  and,  (ii)  a  small  signal  linear  feedback  multivariable  regulator.  The 
equations  of  motion  are  linearized  about  the  trim  condition  of  wings  level,  constant 
altitude  flight,  which  simplifies  the  dynamics  to  allow  separation  between  the  lateral  and 
longitudinal  control  channels.  The  linear  compensator  also  entails  integral  action.  Thus, 
the  controller  consists  of  a  strongly  nonlinear  feed-forward  module  and  a  linear  small- 
signal  multivariable  tracking  controller.  The  novel  proposed  multivariable  set  point 
controller  encompasses  full  three-axes  autopilot  functions.  Moreover,  this  controller  is 
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used  as  a  tracking  controller,  provided  that  the  bandwidth  of  the  command  signal  is 
substantially  less  than  the  bandwidth  of  the  closed-loop  flight  control  system,  and  we 
then  refer  to  the  phugoid  damping  controller.  It  is  envisaged  that  the  latter  will  receive 
inputs  from  a  higher  level  supervisory  control  module,  thus  affording  UCAVs  a  high 
degree  of  control  autonomy. 

The  paper  is  organized  as  follows.  In  Chapter  2,  the  aircraft  model  is  introduced. 
The  aircraft  parameters  are  representative  of  an  F-16  class  aircraft  flying  at  200  m/s  at 
20,000  ft,  and  are  given  in  Table  1.  Chapter  2  also  introduces  the  wind  axes  used  in  the 
design,  the  nonlinear  equations  of  motion,  and  the  nondimensionalization  used  to 
simplify  the  dynamics.  The  required  aircraft  trim  control  settings  for  a  commanded  trim 
state  are  derived  in  Chapter  3.  In  Chapter  4,  the  equations  of  motion  are  linearized  about 
the  trim  state.  These  linearized  equations  of  motion  and,  in  particular,  the  trim  condition 
of  wings  level,  constant  altitude  flight,  with  trim  speed  equal  to  initial  speed,  simplify  the 
dynamics  and  allow  decoupling  between  the  lateral  and  longitudinal  control  channels. 
These  results  are  used  in  Chapter  5  to  design  small  signal  linear  multivariable  regulators 
for  the  two  control  channels.  The  linear  compensators  also  entail  integral  action.  Thus, 
the  controller  consists  of  a  strongly  nonlinear  feed-forward  module  and  a  linear 
multivariable  small-signal  compensator.  Chapter  6  addresses  issues  impacting  the 
operational  envelope  of  the  controller.  The  novel  multivariable  set  point  controller 
encompasses  full  three-axes  autopilot  functions.  Moreover,  this  controller  is  used  as  a 
dynamic  phugoid  damping  controller,  provided  that  the  bandwidth  of  the  command  signal 
is  substantially  less  than  the  bandwidth  of  the  closed  loop  flight  control  system.  Chapter 
7  presents  the  results  of  extensive  simulations  for  the  phugoid  damping  controller’s 
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performance  evaluation.  The  wide  operational  envelope  of  the  controller  is 
demonstrated.  Concluding  remarks  are  made  in  Chapter  8. 
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2.  WIND  AXES  AND  EQUATIONS  OF  MOTION 


The  aircraft’s  point-mass  equations  of  motion  are  derived  using  a  translating  and 
rotating  frame  of  reference  collocated  with  the  instantaneous  position  of  the  aircraft.  The 
rotating  frame  of  reference,  shown  in  Fig.  1,  is  a  triad  of  wind  axes  defined  as  follows: 
The  x-axis  is  aligned  with  the  aircraft  velocity  vector,  the  z-axis  is  aligned  with  the  lift 
vector,  and  the  y-axis  (out  of  the  wing)  completes  the  right-handed  coordinate  system. 
Side  force  F  is  defined  in  the  direction  of  the  y-axis.  The  wind  axes  frame’s  attitude  is 
specified  by  the  v;/,  y,  and  (])  Euler  angles,  corresponding  to  the  aircraft  "heading",  "flight 
path  angle"  (in  the  pitch  plane),  and  "roll  angle".  For  the  sake  of  clarity,  the  final  rotation 
about  the  x-axis  through  the  angle  0  is  not  explicitly  shown  in  Fig.  1.  Also,  velocity 
vector  rolls  are  envisaged,  for  the  roll  maneuver  is  performed  about  the  velocity  vector 
axis. 
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The  three  angles  y,  and  (j)  are  not  unlike  the  Euler  angles  -  they  position  the 
wind  axes  frame  x,  y,  z  relative  to  the  inertial  frame  Xi ,  Yi ,  Zi .  Note,  however,  that 
while  the  heading  angle  \\i  and  the  flight  path  angle  y  are  state  variables,  the  bank  angle  ^ 
is  a  control  variable.  Also,  note  the  polarity  of  the  flight  path  angle  y,  as  indicated  in 
Figure  1. 

2.1.  Aerodynamic  Angles  Definition 


Figure  2:  Sideslip  Angle  Definition 
In  this  paper,  wind  axes,  rather  than  body  axes,  are  used.  However,  a  discussion 
of  body  axes  is  required  to  define  the  aerodynamic  angles  properly.  Thus,  the  body  axes 
triad  (Xb,  Yb,  Zb)  is  related  to  the  wind  axes  (x,  y,  z)  as  follows:  Initially,  the  body  axes 
are  aligned  with  the  wind  axes.  First,  a  rotation  of  P  degrees  about  the  z  wind  axis  is 
performed.  This  is  followed  by  a  rotation  of  a  degrees  about  the  Y^  body  axis.  The 

broken  line  X^  in  Figure  2  is  the  projection  of  the  Xb  axis  onto  the  (x,  y)  plane.  We  note 

that  this  definition  of  aerodynamic  angles  is  not  the  standard  definition  used  when  body 
axes  and  rigid  body  dynamics  are  used  to  describe  the  aircraft’s  motion.  At  the  same 
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time,  the  above-defined  aerodynamic  angles  correspond  to  the  aerodynamic  angles  used 
in  wind  tunnel  work. 

Note  that,  although  the  reference  frame  rotates  as  an  actual  aircraft,  the  aircraft 
model  is  a  point  mass.  Since  a  point  mass  model  is  used,  moments  are  not  included  in  the 
analysis. 

2.2.  Equations  of  Motion 

The  aircraft  dynamic  model  used  in  this  paper  exclusively  employs  the  aircraft’s 
“slow”  states  and  is  therefore  the  “right”  model  for  model  based  outer-loop  controller 
design,  viz.,  altitude-hold,  Mach-hold,  and  heading-hold  autopilot  design,  and  for 
phugoid  damping  controller  design. 

The  point  mass  equations  of  motion  -  see,  e.g..  Fig.  1  and  [5]  -  are: 


V  =  ^^  +  gsinY,  V(0)-Vo 
m 


Lcoscb  +  Fsintl)  g 

Y  = - - 'L  +  ^cosy,  Y(0)  =  Yo 

mV  V 


-Lsm<.^Fcos» 
mV  COSY 


where  the  forces  are: 


_ 

L  =  q— jS^a^oc 
*0 


T  =  mgp 


D  = 
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F  =  qSs,a,P 

'^o 

and  where  the  nominal  dynamic  pressure  is: 

q=ipV“ 

and  Vq  is  the  initial  velocity. 

Nondimensional  variables  and  parameters  are  introduced  as  follows: 


V:= 


_y_ 

Vo 


-  Vo_ 
(0  2.(0 

g 

q:=^ 

m-g 


Cl 


I 

q 


K:=KCl 


a 


w 


a 


a 

d 


P 


a 


^S.  V 


A 


■n 

J 


•P 


Similarly,  the  (barred)  trim  controls  are  scaled  according  to 
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_  a 
a:=  — 
a 


P:=^ 

a 


^  e  Y  „  ^ 


V  A”W  ) 


w 


p 


Remark  For  constant  altitude  and  wings  level  flight,  the  trim  lift  coefficient 
Cl  =  Cl  .  Hence,  the  trim  lift-over-drag  ratio  is  given  by: 


1 


_ 

The  pertinent  aircraft  parameters  are  given  in  Table  1. 


Table  1:  Aircraft  Parameter  Values 


Wing  Lift  Curve  Slope 

aw 

5.3 

/rad 

Tail  Lift  Curve  Slope 

St 

5.3 

/rad 

Wing  Aspect  Ratio 

AR 

3 

- 

Wing  Span 

b 

9.14 

m 

Zero-Lift  Drag  Coefficient 

Cdo 

0.015 

- 

Acceleration  of  Gravity 

g 

9.81 

m/sec^ 

Parabolic  Drag  Polar  Constant 

K 

0.1118 

- 

Mass 

m 

11,336.4 

kg 

Dynamic  Pressure 

q 

12500 

kg/(m-s^) 

Wing  Surface  Area 

Sw 

27.87 

VCi 

Tail  Surface  Area 

St 

5.086 

m^ 

Initial  Velocity 

Vo 

200 

m/sec 

Efficiency  Factor 

h 

0.95 

Atmospheric  Density  (@  20,000’) 

P 

5/8 

kg/m^ 

The  nondimensional  parameters  and  variables  are  given  in  Table  2. 
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Table  2:  Nondimensional  Parameters  and  Variables 


Actual 

Nondimensional 

Cl 

- 

0.319 

eje^ 

- 

12.1 

K 

0.112 

0.0357 

q 

12,500  [kg/(m-s^)] 

3.13 

t 

t  [s] 

0.0491 1 

V 

V  [m/s] 

0.05-V 

a 

a 

16.6- a 

a 

- 

0.0602 

p 

p 

3.03- P 

CO 

(5  [s  '] 

20.4- (5 

0.0694  [s'*] 

4~2 

Using  the  above  parameterization  in  Eqs.  (1),  the  elegant  nondimensional  aircraft 
equations  of  motion  are  derived: 

Nondimensional  Equations  of  Motion 


V  =  sinY-qCDoV'+p-KV"a"  ,  V(0)  =  1 


(2) 


Y  =  V(a cos (t)  +  P sine]))  ,  Y(0)=Yo 


V 


(3) 


V  ^(-asin(l)  +  Pcos(t))  ,  vi/(0)-0  , 


(4) 


COSY 


The  two  parameters  in  the  equations  of  motion  (2)  -  (4)  are  (qC^p)  and  K,  where,  for  the 


flight  condition  considered  in  the  research,  the  numerical  values  are  qC^^  =  0.04695  and 


K  =  0.0357. 

2.3.  Chapter  Summary 

This  chapter  introduced  the  wind  axes  used  to  derive  the  aircraft’s  point-mass 
equations  of  motion.  This  is  a  rotating  frame  of  reference  collocated  with  the  aircraft’s 
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center  of  gravity  (CG).  The  aerodynamic  angles  used  in  the  research  were  defined, 
relating  the  wind  axes  to  the  inertial  and  the  body  axes.  The  aircraft  parameters  and 
variables  were  introduced,  followed  by  the  fundamental  equations  of  motion  and  a  series 
of  nondimensionalizing  variable  and  parameter  definitions.  These  were  used  to  derive 
the  nondimensional  equations  of  motion  used  in  the  rest  of  this  thesis. 
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3.  AIRCRAFT  TRIM  EQUATIONS 

The  equations  of  motion  reported  in  the  previous  chapter  define  velocity  V,  flight 
path  angle  y,  and  heading  \|/  in  terms  of  thrust  setting  |X,  angle  of  attack  a,  bank  angle  (|), 
and  sideslip  angle  p.  In  this  research,  (V,  y,  \|/)  will  be  designated  as  states,  controlled  by 
the  control  variables  (fi,  a,  (|),  P).  In  order  to  use  these  control  variables  to  bring  about  a 
desired  state,  it  is  necessary  to  know  to  what  values  the  control  variables  should  be  set. 
However,  as  the  task  of  back-solving  a  set  of  nonlinear  differential  equations  to  obtain 
inputs  for  completely  arbitrary  outputs  is  a  daunting  one,  some  simplification  is  in  order. 

The  control  system  will  be  based  around  the  idea  that  steady-state,  or  "trim", 
states  are  the  desired  conditions.  This  is  useful  in  making  the  mathematics  tenable 
because  it  forces  the  state  derivatives  to  equal  zero,  changing  the  nonlinear  differential 
equations  into  nonlinear  equations. 

Two  trim  states  will  be  evaluated.  The  first,  "Trim  1",  is  the  simplest.  It  requires 
setting  the  derivatives  of  all  three  state  variables  to  zero.  Thus,  the  aircraft  flies  at 
constant  speed,  at  constant  flight  path  angle,  and  at  constant  heading.  The  second  trim 
state,  "Trim  2",  is  similar,  but  the  heading  is  allowed  to  change  at  a  constant  rate.  This 
allows  the  aircraft  to  come  to  a  trim  state  in  which  it  is  maintaining  a  turn. 

3.1.  Trim  1 

New  values  for  trim  states  V  .  Y .  V  are  commanded,  and  the  corresponding  new 
trim  control  settings  are  a ,  ]I ,  <]) ,  P . 

Setting  the  LHS  of  the  differential  equations  (2)  -  (4)  equal  to  zero  yields  the 
algebraic  trim  equations.  The  latter  are  solved,  yielding  the  trim  control  settings. 
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First,  note  that  eq.  (4)  yields 


T  P 

tan  <[)  = 

a 


and,  therefore,  inserting  eq.  (5)  into  eq.  (3)  yields 


-  1  -  A 

a  =  =:r cosy -COSO , 

y2  * 


Next,  inserting  eq.  (6)  into  eq.  (2)  yields 


p:  =  qCpoV^  -siny  +  ^cos^Y'Cos^  (j) 


Finally,  combining  eq.  (5)  and  (6),  we  calculate 


sin(t)  =  V'  ^ 


cosy 


cos  ([)  = , 


,  1- 

rv^Y 

i 

^cosy^ 

/  _\2 


a  =  . 


cosy 


F  =  qCDoV^+icos^y-siny-KP^V^ 


In  short,  given  that 

•  The  new  states  'F ,  y ,  V  are  specified. 

•  The  control  P  is  chosen  arbitrarily 

^  The  controls  ]I ,  a ,  (|) ,  are  given  in  Eqs.  (7). 
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Special  case  (a):  If  we  set  P  =  0  (no  sideslip),  then 


F2  K  . 


p  =  qCDoV"+^cos  Y-siny 


Note:  When  vp"  =  constant,  there  is  no  aerodynamic  justification  for  P  0.  However,  P 
^  0  might  be  justified  where  LO  (Low  Observables)  considerations  are  important. 

Special  case  (b):  If  we  allow  (p  to  be  independently  controlled,  then 


P  cosy- sine]) 


a  ==r- cosy -cost]) 


^2  ,  K  _ 


p  =  qCco  V  +  ^cos  Y  -  sin  Y  -  K  cos  y  •  sin  (p 


Special  case  (c):  If,  in  (b)  above,  we  set  (|)  =  0  (no  bank),  then 


p=o 


_  _  cos  Y 


It  =  qCooV^  +  ^cos^  Y  -  sin  Y 


This  is  exactly  the  same  result  as  special  case  (a). 
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3.2.  Trim  2 


Now,  to  expand  the  set  of  permissible  trim  conditions  to  include  constant-rate 
turns,  command  new  trim  values  \i/(t)  ,y,W .  Again,  the  corresponding  trim  control 
settings  are  a ,  ]I ,  (j) ,  P . 

We  command  new  values  for  the  trim  state:  -  -o  •  t  ^  Y  ^  V 

Note:  Without  loss  of  generality,  we  now  consider  a  turn  to  starboard,  so  that  the  ensuing 
trim  bank  angle  (j)  is  positive. 

Recall 

- 

(0  :=  — (D 

g 

t:=^t 

Vo 

Inserting  these  nondimensionalizing  definitions  into  the  previous  trim  state  definition 
produces  the  new  trim  state  definitions 
\1/  =  -  (0  •  t 

where  both  t  and  6J  are  nondimensional.  The  corresponding  trim  control  settings  are 
]I,  a,  p 

Here  (o  >  0,  viz.,  without  loss  of  generality,  we  consider  a  turn  to  starboard,  so  that  the 
ensuing  trim  bank  angle  ([)  >  0. 

Now,  setting  the  LHS  of  the  differential  equations  (3)  and  (4)  equal  to  0  yields 


16 


(10) 


a  sin  (j)  -  P  cos  (j)  =  =cos  y 

a  cos  ([)  +  P  sin  (])  =^cosy 

This  is  a  system  of  two  nonlinear  equations  in  the  unknowns  a  and  ^ ;  in  Eqs.  (10),  V , 
P  and  Y  are  parameters. 

We  now  embark  on  the  solution  of  the  nonlinear  system  of  equations  (10).  We 
momentarily  consider  the  linear  system  in  the  “unknowns”  sin  (|)  and  cos  (|) .  Thus,  we 
have 


a 

-p' 

'^sin<l)^ 

cosy 

.p 

a 

^COS(|)^ 

‘  J 

The  solution  of  this  linear  system  in  sin  (j)  and  cos  ^  is 

-  1  cosy/— 

Now, 

sin^  (|)  +  cos^  (])  =  1 


1  cos  Y^Q^2^2Y2_|_p2_(_^2^p2^2Y2j_j 


(a^  +  p^f 


2  rr 


^2^p2^£OSjr(l^_2v2) 

^  V  ^ 
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/cos^y 

V 

V" 

( 

1 

— 

cosy 

l  +  co"V" 

V 

1 

r  ~  _ 

cosy 

l+co'V^ 

,2r7 


A 


(ovj^^(l  +  co"V')-P^+P 


-®Vp 


ai) 


Finally, 


|i  =  qCi3oV^  -siny  +  ^cos^  y(i+  a)^V^)-KP^V^ 


Special  case:  P  =  0.  Then 


V' 


sin  (j)  = 


cos<t)  = 


(OV 


Vl+to"V' 

1 

Vl+M^V' 


fi2) 


p  =  qCoo  -  sin  y  +  ^cos^  y(l  +  (O^  ) 


In  the  special  case  w  =  0,  P  ^  0,  the  TRIMl  results  in  Eqs.  (7)  are  recovered.  Indeed, 
TRIMI  (Eqs.  (7))  is  a  special  case  of  TRIM2  (Eqs.  (10)). 

3.3.  Trim  Equation  Summary; 

Given  the  new  trim  conditions  V  ,  y ,  -  q  ,  the  corresponding  trim  control  settings  are: 
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P  =  qV'Cj,o-sinY+^(l  +  (o"V') 


^  li  .  T-r2'tr2\ _ 2  t:  ■Lr\T2.v2  -r^l 


cos^Y-KV^K^^co^ 


l+co"V' 


V^  1  1+(1)V  0;^,  T>r2  tt2 


l+(0  V  COSY  V  ^ 


-cos^Y-K^CO^  +Ke.^«^V 


It  is  remarkable  that  the  aerodynamic  controls  a  and  ^  are  not  dependent  on  the 
problem  parameters  (qC^^ )  and  K.  The  throttle  setting,  jl ,  is  however,  dependent  on  the 

problem  parameters  (qCp_^ )  and  K. 

A  reasonable  choice  for  K^-  is 
K^-  =  0.05 

Indeed,  the  following  holds: 

3.03P  =  Kp_-20.4-(O 

where  P  and  to  are  the  physical  (original)  variables. 

For  (0  =  15° /sec,  use  P  =  5° 

3.03-5  =  Kp- •20.4-15 


3  03 

^  V  ^=^^  =  0.05 
20.4-3 
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Remarks 


•  For  coordinated  "bank-to-tum"  turns  (BTT),  there  is  no  sideslip  angle;  choose 

•  Conversely,  for  skid-to-tum  (STT)  control  there  is  no  bank;  set 

^  =  0 


which  leads  to 


V' 


a  =  T=;rCOS  Y 

^  f  2  * 


o-  0)  _ 

B  =  -=cosY 
V 


(14) 


K 


F  =  CooqV^  -  sin  Y  +  ^cos'  Y 


If  (t)  is  independently  controlled,  e.g.,  in  order  to  point  the  air  vehicle’s  weapons 
for  target  tracking  and  fire  control,  the  trim  equations  are 


cos  ([)  sin  ([) 
-  sin  ([)  cos  (j) 


o'!  1  J  I 


=  ^cosv| 


V 


whereupon  the  solution  is  obtained 
a  =  =^(cos  (])  +  CO V  sin  (l))cos  y 


V 


P  =:^(sin(t)-(oVcos(t))cosY 
p  =  Cpop -  sin  Y  +  ^  (w  V  sin  (j)  +  cos  cos^  Y 


(15) 
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3.4.  Chapter  Summary; 

Given  the  equations  of  motion,  the  required  control  inputs  were  solved  for  a  given 
desired  steady-state,  or  “trim”,  set  of  states.  Two  trim  configurations  were  investigated: 
“Trim  1”,  a  simple  straight-line  flight  path  ( tj/ ,  y ,  V ),  and  “Trim  2”,  a  trim  condition 
that  allowed  steady  turns  (\|/  =  wt ,  y ,  \  ).  The  fact  that  there  are  four  control  variables 
and  only  three  states  means  that  the  system  is  over-determined;  ^  and  P  are  redundant. 
Cases  where  each  was  used  were  examined.  Table  3  shows  the  equation  numbers  for  the 
various  trim  control  laws  presented  herein. 


Table  3;  Trim  Control  Laws  Summary 


P  Controlled 
Independently 

p=o 

P  =  K,.(-(S) 

4)  =  0 

(j)  Controlled 
Independently 

Triml 

V  Y  ¥ 

Eqs.  7 

Eqs.  8 

N/A 

Eqs.  8 

Eqs.  9 

Eqs.  11 

Eqs.  12 

Eqs.  13 

Eqs.  14 

Eqs.  15 
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4.  LINEARIZATION  ABOUT  THE  NEW  TRIM  ^  ,Y,  -  w 

Modem  control  theory,  with  its  vast  assortment  of  mathematical  tools,  is 
predominantly  concerned  with  linear  models.  In  order  to  take  advantage  of  this,  the 
nonlinear  aircraft  model  must  be  linearized.  Linearization  takes  place  about  a  specified 
operating  point.  Two  such  operating  points  are  introduced:  one  straight  and  level,  the 
other,  a  wings-level  climb.  Flight  conditions  entail  a  constant  velocity.  This  choice  of 
trim  points  for  linearization  allows  decoupling  of  the  lateral  and  longitudinal  channels. 

The  differential  equations  in  the  perturbation  states  v,  y,  \\i  and  in  the  controls’ 
perturbations  are: 

v  =  cosY-Y-2qCD(,V-v  +  |j,-2KVa^  •v-2KV^a-a 


Y  =  — V  — (acostj)  —  Psintj))- V  — v(cos(|)-a— asin(t)-(l)  +  sin(t)'P  +  Pcos(|)-(|)) 


xj/  =  — ^^(-asin(t)  +  Pcos(t))-  v  +  -^^^^^^(-asin(|)  +  Pcos(t))-Y 
cosy  cos  Y 


+  - 


cosy 


—  (-  sin  (])  •  a  -  a  cos  (])  •  (|)  +  cos  (j)  •  P  -  P  sin  (j)  •  (])) 


Use  trim  equation 


_  _  _  COS  Y 

a  cos  (|)  +  P  sin  (])  =  _ ^ 


Y  =  .  Y  -  2^2^  •  v  -  V  cos  (|)  •  a  -  V  sin  (j)  •  P  +  v(a  sin  (j)  -  P  cos  (j))-  (j) 


V 


V 


Now,  trim  relations  also  yield 


—  •  A  Q-  A  wcosy 
a  sm  ([)  -  p  cos  0  =  — = — 


COS  Y  sin  Y  _  —  _  —  —  —  n 

Y  =  -2  — ^  •  V  — ^  ’Y-Vcostlx-a  +  O)  cos  y  •  ([)  -  V  sin  (j)  •  p 


22 


W  --  V.;r  l.V 

\1/  =  -=•  v-cotanv-Y - -sin  (p-a-=-(|)  +  — ::Cos(])-p 

V  cos  Y  V  cos  Y 

Hence,  the  linear  perturbations  equations  are 


dt 


V 

V 

a 

Y 

=  A 

Y 

+  B 

y. 

y. 

.P. 

where  the  dynamics  matrix  is: 

COSY  0 

sin?  Q 
V 

(otanY  0 


0 

-  V  sin  (|) 

yCOS^ 

COSY  _ 

Using  the  expression  for  the  trim  AOA  (angle  of  attack)  a  ,  we  calculate  the  A  and  B 
matrices’  parameters 

qC,„  +  Ka“  =qC„  +K(l  +  ffl"V")^^-Kp 

and 

KV^cc  =  K-Jcos^  Y  •  (l  +  cb^V" )- 

Similarly,  inserting  the  expression  for  $  ,  we  calculate  the  additional  B  matrix  parameters 


-2v(qCi,o+_Ka") 

qCosy 

_  0) 

~V 


and  the  control  matrix  is: 


B  = 


-2KV^a 

-  V  cos  ([) 
sin  ({) 


-V- 


cosY 


0 

mcosY 

V 
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-  -  V'  1 

V  sin  (j)  =  — —5=7 - z 

1  +  0)  V  cosy 


0)vj^^^.^  cos^  Y  -  +  P 


J 


V  cos  (j)  : 


1 


^  ll  +  co^V 


1+  o)^V^  cosy 


cos^y-P^  -o)VP 


4.1.  Linearization  About  Wings  Level  Climbing  Flight,  Trim  Speed  =  Initial  Speed 

Now,  for  the  special  case 

V=l,  ^  =  0,  P=0,  co  =  0 
we  calculate  the  parameters 

qCpo  +  Ka'  =  qCco  +  Kcos"  y 


KV^a  =  Kcosy 


Hence,  the  dynamics  matrix 


A  = 


■2(qC^o 
-2cosy 


+  Kcos^  y) 


cos  Y  i  0 
-siny  I  0 
0 


-2- 


cosy  I  0 


-2cosy  -siny  i  0 

"  ""o . . . 6  ”  To 


and  the  control  matrix 


1 

-2Kcosy  ! 

0 

0 

B  = 

0.. 

0 

0 

0 

0  i 

1 

-1 

1 

cosy 
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Note  that  the  block  diagonal  structure  of  the  A  and  B  matrices  renders  the  longitudinal 


channel  decoupled  from  the  lateral  directional  channel  \|/ . 

UJ 

4.2.  Linearization  About  Trimmed  Win2s  Level.  Constant  Altitude  Flight,  Trim 


Speed  =  Initial  Speed 


Now,  for  the  important  special  case 

V=l,0  =  O,Y  =  O,P=O,a  =  l,(o  =  O 
we  obtain  the  dynamics  matrix 


a 

O 

+ 

1 

0 

1 

N) 

o  I 

O  I 

0 

o 

O 

and  the  control  matrix 


'1 

-2K 

0 

B  = 

0 

-1 

0 

0 

0 

-1 

Moreover,  it  is  readily  verifiable  that 
qC„„  +  K  =  Ss- 


where  and  Cl  are  the  lift  and  drag  coefficients  for  wings  level,  constant  altitude, 
trimmed  flight. 

Hence,  the  dynamics  matrix  is 

■-2/(q,/C„)  liO" 

A=  -2  0  i  0 
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Thus,  the  phugoid  dynamics  are  exclusively  determined  by  the  solution  of  the 
characteristic  equation 


i.e.. 


Discussion: 

The  nondimensional  natural  frequency  is 

(0n=V2, 

the  dimensional  natural  frequency  is 


and  the  damping  ratio  is 


We  see  that  the  higher  the  lift  over  drag  ratio,  the  lower  the  damping  of  the  phugoid  is. 
Moreover,  the  phugoid  is  oscillatory.  The  phugoid  ceases  to  be  oscillatory  for  very  low 


lift  to  drag  ratios 


Cl  1 

=^<-r=- 

Cd 
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—  —  c 

Note:  In  our  case,  Cl  =0.319  ,  Cj,  =  0.0264  ,  and  =^  =  12.1 
4.3»  Chapter  Summary 

In  order  to  design  the  (model-based)  small  signal  controller,  the  nonlinear  aircraft 
model  was  linearized.  Linearization  takes  place  about  a  specified  operating  (trim)  point. 
Two  such  operating  points  were  introduced:  one  straight,  level,  constant  velocity,  the 
other,  a  level  climb  with  constant  velocity.  This  choice  of  linearization  allowed 
decoupling  of  the  lateral  and  longitudinal  channels. 
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5.  LINEAR  CONTROL  DESIGN 


The  aircraft  model  was  linearized  in  the  previous  chapter.  Thus,  a  linear  small- 
signal  controller  can  be  designed.  The  small-signal  linear  controller  is  driven  by  error 
signals.  Since  the  nonlinear  plant  will  be  operated  about  a  new  set  point,  a  nonlinear 
feed-forward  controller  is  used  in  parallel  with  the  linear  feedback  portion.  The  nonlinear 
feed-forward  module,  the  “trim  solver”,  is  based  on  the  aircraft  trim  equations  derived  in 
Chapter  3.  Proportional  controllers,  using  the  linearized  plant  dynamics,  are  designed 
separately  for  both  lateral  and  longitudinal  channel  control,  and  integral  action  is 
included  in  the  lateral  channel.  Three  separate  equivalent  gain  options  are  calculated  for 
the  longitudinal  channel’s  small-signal  controller,  referred  to  as  options  (la),  (lb),  and 
(2). 

5.1.  Theory 

Consider  the  nonlinear  plant 

x  =  f(x,u),  x(o)  =  Xo,  x,Xog91",  ugSI™,  0<t 
y  =  x 

and  the  reference  signal 
r  =  X 

where  x  is  a  rest  point  (or  equilibrium  point,  or  trim  point)  of  the  dynamical  system.  In 
other  words,  3  u€  s.t.  f(x,u)  =  0 . 

Consider  the  set  of  trim  points  of  interest  X  c  91" . 

Hypothesis: 

V  xeX3uG9T"  s.t.  f(x,u)  =  0. 
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Moreover,  under  some  mild  assumptions  on  the  function  f(x,u),  the  implicit  function 
theorem  yields: 

3  smooth  function  g  :  91"  — >9^'"  s.t.  u=  g(x) 

In  other  words, 
f(x,g(x))  =  0  VxgX 

i.e., 

f  (x,  g(x))  =  0  { end  of  hypothesis) 

Given  the  reference  signal  r  =  x  ( =  constant),  consider  the  control  law 
u  =  g(r)  -  V 

where  v  is  the  output  of  the  linear  controller.  The  dynamics  then  are 
x  =  f(x,g(x)-v)  x(o)  =  Xo  0<t 
y  =  x 

WeTl  linearize  the  function  f  about  the  trim  point  (x,g(x)).  Thus, 
f  (x,  g(x)  -  v)  =  f  (x,  g(x)) + A(x  -  x)  -  Bv  +  H.O.T. 

=  0  + A(x-x)-Bv  +  d 

where 

A  =  —  ,  B  =  —  ,  d  =  H.O.T.  (higher  order  terms) 

(s,g(x))  (x.g(x)) 

Hence 

-e=-Ae-Bv-d 

where 

e=r-y =x-x 
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Therefore,  the  perturbation  dynamics  are 
e  =  Ae  +  Bv  +  d 

Tracking  of  a  constant  command  signal  r  is  guaranteed,  provided  that  the  control 
signal  V  is  generated  by  a  stabilizing  state  feedback  control  law.  We’ll  use  a  PI  control 
law.  We  need  integral  action  because  we  need  to  reject  the  linearization-induced 
disturbance  d  [9].  Moreover,  although  we  do  not  need  integral  action  for  tracking  the 
constant  set  point  command  when  the  linearized  dynamics  with  d  =  0  are  used,  integral 
action  might  also  help  with  tracking  a  dynamic  reference  signal. 

Hence,  we  augment  the  perturbation  dynamics  as  follows; 
e  =  Ae-l-Bv  +  d 


z  =  e 


where  z  is  the  “charge”  on  the  integrator. 
Use  the  linear  control  law 

v  =  Kpe  +  KjZ 

In  matrix  form,  closed-loop  dynamics  are 


e 

1 — 

o 

< 

e 

'b' 

r 

+ 

V  + 

z 

I  0 

z 

0 

0 

and  the  linear  controller  output  is 


v  =  [K,  K,] 


By  substitution. 


d 

A  +  BKp 

BK,' 

r 

zn 

+ 

dt 

z 

I 

0 

z 

0 

^  ) 

V  J 

e(o)  =  eo 
z(0)  =  0 
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^  A  +  BKp  a  stability  matrix.  Again,  the  poles  of  the  pair  (A,  B)  are  assignable 
because  (A,  B)  is  a  controllable  pair  (in  view  of  the  fact  that  the  nxn  matrix  B  is 


nonsingular.) 

Hence,  we  first  design  the  proportional  controller.  This  “governs”  the  tracking 
error  dynamics.  The  second  step  entails  the  setting  of  the  integral  gains. 
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Note:  Integral  action  is  not  strictly  needed  when 

a)  There  are  no  disturbances,  viz.,  the  linearized  dynamics  are  used,  and 

b)  The  reference  signal  r  =  x  =  constant. 


The  phugoid  damping  control  system’s  block  diagram  is  shown  in 

Figure  3. 


Nonlinear  Controller 


Figure  3:  Phugoid  Damping  Control  System 

Here 

y  =  x 
e  =  r-  y 
z  =  e  ,  z(0)  =  0 

u  =  g(r)-Kpe-KiZ 
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where,  we  recall, 


fv^ 

a 

Y 

,  x  = 

Y 

,  u  = 

<1) 

w 

V 

P 

Zooming  in  on  the  \|/-channel  -  see  Figure  4: 


Figure  4:  v|i-Channel  Block  Diagram 

The  function  g(»)  implements  the  nonlinear  trim  signal  calculation  developed  in  Chapter 
3.  The  nonlinear  trim  function  g(»)  is  explicitly  specified  as  follows: 


a 


:g(v,Y,(i);P) 


where  (see  Eqs.  (11 )) 


P  =  qCcoV"  -  sin  Y  +  ^cos"  y(i  + 
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•  A  ^ 

sin(|)  = 


1 


cosy  1+co^V^ 


(0V^^^(l  +  (0^V^)-P^+3 


COS(t) 


cosy  1+co^V^ 


,2  T7 


^^(l+co^V^)-|3^-(oVP 


5.3.  Proportional  Control 

The  controller’s  design  is  based  on  the  linearized  special  case  of  trimmed,  wings 
level,  constant  altitude  flight,  with  trim  speed  =  initial  speed.  Control  of  the  lateral  and 


longitudinal  channels  is  addressed. 
5.3.1.  Lateral  /  Directional  Channel 


The  desired  flight  control  system  (PCS)  lateral  channel’s  aircraft  control 
parameters  are  specified  in  Table  4. 

Table  4:  Desired  Parameters  of  Lateral  PCS  Channel 


Dimensional 

Nondimensional  j 

T 

10  [sec] 

0.49 

-Jljl 

V2/2 

“n. 

^0 

Note:  The  lateral/directional  control  channels’  nondimensional  time  constant 


In  the  lateral  /  directional  channel,  the  reference  signal  is  \|/ ,  and  thus  the  error  signal  is 
e^=\j/-\l/  [rad] 
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Now,  the  linearized  lateral  directional  equation  of  motion  is 


S  =  -(t>  +  P 

The  controls  are 


4*  ~ 


P  =  K,,, 


e 


¥ 


^P,v) 


e^(t)  =  e 


e,,(0) 


We  want 


e,;,(t)  =  e  "e,^(0) 

Hence,  set 

Note  that  x  is  the  specified  nondimensional  directional  time  constant. 

Choosing  a  dimensional  x  =  10  sec  ^  nondimensional  x  =  0.5. 

K,.^-Kp,,,=2  (16) 

If  a  BTT  control  strategy  is  employed,  then  Kp  =00 

K,,,,=2  (17) 

If,  instead,  an  STT  control  strategy  is  employed,  then  =  0 
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■2 


(18) 


^  Kp,^=- 


In  conclusion: 

1 .  The  design  specification  is  the  directional  channel’s  time  constant  t. 

2.  We  operate  on  e^ 

3.  The  control  variables  are  (j)  and/or  p. 

4.  The  actual  control  signals  are: 


(t)  =  -K^,^e^  +  Asin 


1 


cosy  1+co^V^ 


,2r. 


=  -2e^  +  Asin 


cosy  l  +  co^V^ 

P  =  -Kp,^e^-Kp-CO 


2  —  . 

coVj 

V 

yA  \  / 

r_(0 


yj 


(19) 


5.  In  pure  BTT  control,  Kp^  =  Kp  _  =  0  ,  viz.,  P  -  0. 

6.  In  pure  STT  control,  (])  =  0. 

5.3.2.  Longitudinal  Channel 

The  reference  signals  are  V  and  y ,  and  thus  the  error  signals  are 


ev=V-V 

=  y-y  [rad] 

which  are  fed  back  to  the  linear  controller.  Thus,  the  linear  controller  operates  on 
feedback  signals  Cv  and  Cy .  The  control  variables  for  the  aircraft  longitudinal  channel  are 
a  and  p  .  We  need  to  specify  the  four  controller  gains  Kav,  Kay,  K^v,  K^Y  .  We  will  use 
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(longitudinal  channel’s  natural  frequency)  and  (longitudinal  channel’s  damping 

ratio)  as  the  longitudinal  channel’s  design  specs. 

The  design  equations  are  derived  as  follows: 

The  closed-loop  system’s  dynamics  matrix  is 

-2S  +  K,,-2K-K„  I  +  K„-2K-K, 

‘'-2-K.v  -K^ 

The  characteristic  equation  for  a  2x2  matrix  is 

0  =  A.'+2^p(0„^+0)^^ 

=  }^-  trace(AcL^ )  +  det(AcL_, ) 

Hence,  the  design  equations  are 

Trace(  Acl^  )  =  -2^pa)„^ 

Det(  Aqlp  )  — 

Thus,  the  design  equations  yield  two  equations  in  four  unknowns  (the  gains  Kav,  Kay, 

K|dV5  K|j,y)- 

( 

Kccy  +  2K  •  K^y  —  Kj^y  =  2  “ 

V 

K„f  2K  ■  K„v  -  K„  +  2^1+  (2  +  K„  Xl  +  K„  -  2K  ■  K„  )=  »;_ 

I  '-L  J 

The  first  equation  yields 

2K  •  +  2^  =  24p(0„^  - 

Inserting  this  expression  in  the  second  equation  yields 
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KaY(24p<Mnp  •■K„J+(2  +  K„v)(i  +  K^^  -2K-K„J-( 


Hence,  we  have  obtained  the  two  design  equations 


K  +2KK^y  K  y  2  — 


-ay  ‘  ^  'oP'^np  ^ 


K„(24,<0.;  -  )+  (2  +  K„  Xl  +  K„  -  2KK^  )=  <2») 

In  summary,  we  have  two  equations  in  the  four  unknowns  Kav,  Kay,  K^v,  K^y  •  In 
the  design  equations,  only  the  problem  parameter  K  features  —  and,  of  course,  the  design 
specs  and  . 

Finally,  the  complete  control  law  is 


p,  =  qV^Cpo  -  sin  Y  +  i  (l  +  (O^V^  )cos^  Y  “  KV^Kp-(0^  -  •  e^  -  •  e^ 


(l  +  (D  V  )  Kp  _co  K„y -ey  Kjjy -e^ 


I  1+(0^V^  2-  x.2  -2  -  ir 

<^  =  Asm]  - r  wVJ — =4 — cos  Y-Kp^gO)  -%  co 

1+co  V  COSY  V  V 


P  ^6.0)  ( 


Obviously,  when  Kav  =  Kay  =  K^v  =  K^y  =  0  ,  we  are  back  to  the  open  loop  case.  Hence, 
the  design  equations  yield 


ty„  =V2 

rip 
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=  V2  =  a)„ ,  and  ^ ,  as  expected. 

’’  V2  Cl 

1.  Choose  Kn,v  =  0  .  K,,y  =  0 

This  set  of  gain  restrictions  is  somewhat  arbitrary.  It  is  chosen  because  it  is 
intuitively  obvious.  Within  the  controller,  velocity  error  has  no  effect  on  AO  A,  and  flight 
path  angle  error  has  no  effect  on  thrust  setting.  This  implies  that  AOA  is  affected  only  by 
an  error  in  flight  path  angle  Cy,  while  thrust  setting  is  controlled  solely  using  velocity 
error  ey.  This  is  similar  to  the  logic  that  a  pilot  might  use  to  maneuver  an  aircraft  in 
certain  situations,  using  pitch  to  control  the  direction  of  the  aircraft  and  the  throttle  to 
control  the  speed.  Thus,  the  design  equations  are: 


Kay-K,,=i^p(0„^-^ 


2^,a)„K^  -Kl^+2- 4KK^  =  co^ 

Thus,  we  have  two  equations  in  the  two  unknowns  Kay  and  K^v  .  The  second  design 
equation  is  a  quadratic  equation  in  Kay : 

K^,+2(2K-^p(0„>„,+(o;  -2  =  0 


Kay  =  -  2K  ±  +  2  +  4K^  -  4K^p(0„ 


K,,  =  2^  -  2K  -  ^p0)„^  ±  M  -  iKp  +  2  +  4K^  - 


If  we  specify 
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then  the  two  solutions  are 


We  want  to  increase  the  bandwidth  of  the  closed-loop  flight  control  system  in  order  to 
enhance  tracking  performance.  Hence,  we  will  specify 


0)„,>0)„=V2 

Remark  1 :  If  we  should  choose  =  &)„  =  V2 ,  then  in  case  (a) 
=2(1- 2K)  =  1.8572  =  2 

fc  ^ 

K,v=2^-2K  =0.0226=0 

J 

and  in  case  (b) 
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This  choice  of  cOn  and  ^  effectively  changes  the  damping  of  the  closed-loop  system  while 
retaining  the  natural  frequency  found  in  the  open-loop  system.  Since  settling  time  is 
related  to  natural  frequency,  system  response  can  be  improved  by  increasing  the  natural 
frequency.  Therefore,  we  will  arbitrarily  increase  the  desired  natural  frequency  to 
(0^^  =  Vs  .  Note  that  could  be  increased  further,  at  the  expense  of  increasing  the 
gains  of  the  controller  matrix. 

Note  also  that,  in  each  case,  three  of  the  four  gains  are  set  to  approximately  zero. 
These  are  particularly  low-gain  controllers.  This  is  not  surprising,  given  the  limited 
design  requirements  given. 

Remark  2:  A  physically  realizable  controller  matrix  is  required  to  be  composed  of  all 
real  parts.  A  negative  discriminant  would  produce  complex  gains.  We  need  a 
nonnegative  discriminant,  viz.. 


— (0„  -2k1  >0)^  -2 

2  “p 


0<(O„  <2V1  +  4K 


■V2k) 


Hence,  choose,  e.g.. 


a)  K„  =4=fVl-4V6K  +  8K'  -tV3-2V2K 


K,v  =2&  +  ^^Vi-4V6K  +  8K'  -V3-2V2K 
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Thus,  Kcey  >  0  and  K^v  <  0,  as  required. 


Again,  Kay  >  0  and  K^v  <  0,  as  required.  Now  the  gain  Kay  is  lower  than  before,  and  the 
gain  K^v  is  higher. 

2.  Choose  K»y  =  0  ,  K|iv  =  0 

Now,  within  the  controller,  velocity  error  has  no  effect  on  thrust  setting,  and  flight 
path  angle  error  has  no  effect  on  AO  A.  This  implies  that  AOA  is  affected  only  by  an 
error  in  velocity  error  ey,  while  thrust  setting  is  controlled  solely  using  flight  path  angle 
ey.  This  is  the  opposite  of  the  previous  option,  and  may  also  be  used  naturally  by  pilots 
flying  in  certain  situations.  Then  the  gains  are: 
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If,  in  addition,  we  specified  a)„  =  ^^2,  to  maintain  the  open-loop  natural  frequency  in  the 


closed-loop  system  while  increasing  damping  ratio  to  the  specified  ,  then 


K 


^-1 
2K  +  1-S^ 

Cl 


Choosing,  as  before, 

“np 


Concerning  the  polarity  of  the  gains:  We  see  that  when  an  increase  in  airspeed  is 
required  (ev  >  0),  the  control  law  brings  about  a  reduction  in  the  angle  of  attack,  as 
expected.  However,  when  an  increase  in  flight  path  angle  (toward  a  steeper  dive)  is 
required,  the  control  law  brings  about  an  increase  in  thrust  setting,  i.e.,  the  polarity  of 
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does  not  make  physical  sense,  although  the  linearized  closed-loop  flight  control  system  is 
stable;  we  would  prefer  K^y  >  0. 


Demanding  K^y  >  0  requires 


Numerically,  this  is  equivalent  to 
0<o)„  <0.01596 

=  0.0112V2 

which  entails  a  very  low  bandwidth. 

Alternatively, 

(D„  >19.794 

Bp 

-I4V2 

yielding  a  very  high  bandwidth. 

Thus,  given  the  constraints  of  setting  gains  Kay  =  K^v  =  0,  the  intuitively  expected 
gain  sign  configuration  can  be  achieved,  but  only  if  ©n  is  set  to  a  very  high  or  very  low 

level.  Setting  it  to  such  a  low  level  as  0.01 12V2  is  tantamount  to  removing  the 
controller  completely;  phugoid  oscillations  would  continue  for  a  long  time.  On  the  other 
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hand,  setting  a)„  =  Uyf2  could  produce  desirable  results.  Settling  time  would  be 
decreased  greatly.  However,  the  gains  required  for  such  performance  would  be  much 
higher  than  those  calculated  previously.  These  higher  gains  force  the  control  actuators  to 
work  harder,  possibly  to  the  undesirable  point  of  saturation.  Even  worse,  the  system 
response  could  be  increased  to  a  point  where  system  phase  margins  decrease  to  levels 
where  actuator  lag  affects  system  stability. 

It  is  apparent  that  there  is  nothing  “wrong”  with  the  analysis;  all  specifications 
were  met,  and  the  system  is  stable.  These  are  the  gains  required  to  bring  about  the  stated 
^  and  (On  ,  given  that  two  variables  (or  degrees  of  freedom)  were  expended  in  choosing 
Kocy  =  Knv  =  0.  On  the  other  hand,  this  “curiosity”  in  the  gain  polarity  is  reflected  in  the 
system’s  performance.  For  example,  given  a  excess  velocity  error,  the  controller  will 
command  an  increase  in  AOA.  The  gain  involved  in  this  step  is  much  higher  than  the 
others  calculated,  so  the  AOA  increase  is  relatively  high.  The  aircraft  will,  naturally, 
begin  to  climb.  The  curious  gain  polarity  then  will  cause  the  thrust  setting  to  increase. 
This  thrust  increase  further  perturbs  the  dynamics,  exacerbating  the  error.  However,  the 
AOA  increase  is  sufficiently  high  to  force  a  velocity  decrease  in  spite  of  the  thrust 
increase.  The  gains  work  against  each  other,  balanced  in  such  a  way  as  to  produce  the 
desired  ^  and  cOn.  This  is  obviously  inefficient  in  terms  of  actuator  usage,  but  it  is 
unavoidable  unless  the  design  specifications  are  changed. 

5.4.  Inteeral  Action  Control 
Let 


45 


dfe"|_rA  +  BKp  BKjYe^  rijY  e(o)-eo 
dt[zj  [  I3  O3  j^zj  [03J  ’  z(0)=0 
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The  system  is  composed  of  two  independent  channels:  lateral  and  longitudinal.  The  two 
channels  can  be  decoupled  in  the  same  manner  as  was  done  in  the  proportional  controller. 
The  lateral  /  directional  channel  dynamics  are 

Ve  =  -K,,^Ve  -  (0)  =  r^ 

Zv(o)  =  0 

i.e., 

1 

¥.(0)  =  V.. 

2,(0)  =  0 

5.4.1.  Lateral  /  Directional  Channel 
Recall  the  form  of  the  closed-loop  system  description 
e  =  Ae  +  Bv 

[A  +  BKp  BKjl  ^ 

where  A  is  the  closed-loop  linear  dynamics  matrix  =  ^  ^  .  For  this 


Thus,  the  closed-loop  linear  dynamics  matrix  is 
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A  +  BKp  BK, 
I  0 


~  ^  +  Kp_^ 


and  the  characteristic  equation  is 


X"+(Kj„-Kp„)>.  +  K,,.-K|,,.=0 


We  want 


<  =K^.z-Kp,, 


^P.'t 


Choosing,  as  before, 


—JC  = 


-—i= - 

^  42-x 

Choosing  the  same  time  constant  as  was  selected  in  the  previous  proportional  case  would 
allow  the  proportional  and  the  integral  parts  of  the  controller  to  track  similar  dynamic 
reference  signals.  Therefore,  choosing  x  =  14  ,  as  before,  yields 

Hence,  the  integral  action  gains  satisfy  the  equation 

K,,.-Ks,,=2  (25) 

If  a  BTT  control  strategy  is  employed,  implying  Kp  ^  =  Kp  ^  =  0 ,  this  defines  , : 
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K|i,,=Kp,=0 


K„  =  K,,.  =  2 


} 


(26) 


If  a  STT  control  strategy  is  employed,  implying  ,,,  =  ^  =  0 ,  this  defines  Kp  ^ : 


Kp,^  =  Kp,,  =  -2 


} 


(27) 


If  a  hybrid  BTT/STT  control  strategy  is  used,  we  could  somewhat  arbitrarily  choose,  e.g.. 


Kp.^=Kp,,=-0.5 


to  weight  bank  angle  twice  as  heavily  as  sideslip  angle  in  the  controller. 

5.3.2.  Summary 

When  integral  aetion  is  used  in  the  lateral  /  directional  channel,  the  complete  control  law 
is  augmented  as  follows: 


([)  =  A  sin 


1 


1+  (o^V^  cosy 


|l+(O^V\__2_  -2 


mV, 


cos^Y-Ki^w^-K^^O) 


,0) 


P~  Kp  _-(0  Kp^e^  Kp^^ 

5.5.  Chapter  Summary 

Since  the  aircraft  model  was  linearized  in  the  previous  chapter,  a  linear  controller 
could  be  implemented  in  this  chapter.  The  theory  of  the  control  system  was  introduced, 
with  the  controller  driven  by  error  signals.  Since  the  plant  was  nonlinear  and  had  to  be 
operated  about  a  nonzero  set  point,  a  nonlinear  feed-forward  design  was  implemented  in 
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parallel  with  the  linear  feedback  portion.  This  feed-forward  portion,  the  “trim  solver”, 
was  based  on  the  aircraft  trim  equations  derived  in  Chapter  3.  Proportional  controllers 
were  designed  separately  for  both  lateral  and  longitudinal  channel  control,  and  integral 
action  is  added  to  the  lateral  channel.  Three  separate  equivalent  gain  options  were 
calculated  for  the  longitudinal  channel,  referred  to  as  options  (la),  (lb),  and  (2).  A 
summary  of  the  gains  used  in  the  linear  compensator  is  given  in  Table  5.  The  nonlinear 
portion  of  the  control  law  was  summarized  previously  in  Table  3. 


Table  5:  Linear  Control  Gains  Summary 


\ 

Lateral 

Longitudinal 

Generic 

BTT 

STT 

Koty~K|xv— 0 

(Option  la) 

Koty— K^V— 0 

(Option  lb) 

KaV— Kfiy— 0 
(Option  2) 

Proportional 

Eq.  16 

Eq.  17 

Eq.  18 

Eqs.  23 

Eqs.  24 

Eqs.  25 

Integral 

Eq.  26 

Eqs.  27 

Eqs.  28 

N/A 

N/A 

N/A 
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6.  DYNAMICS  DEPENDENCE  ON  OPERATING  POINT 


The  phugoid  dynamics  aircraft  model  under  consideration  is  strongly  nonlinear. 
Since  high  amplitude  maneuvers  and,  conversely,  large  excursions  in  the  state  variables 
are  contemplated,  it  is  most  instructive  to  investigate  the  linearized  dynamics’ 
dependence  on  the  state  variables. 

6.1.  Drae  Polar 

At  trim,  the  aircraft’s  “drag  polar”,  or  plot  of  drag  vs.  velocity,  is  equal  to  the 
thrust  required  to  maintain  that  trim  state.  This  thrust,  calculated  previously,  was  found 

to  be  jl  =  CooqV^  -  sin  y  +  ^  (w  V  sin  ^  +  cos  cos^  y  •  At  the  trim  conditions, 

Y  =  w  ^  =  0 ,  thrust  p  varies  with  velocity  as  shown  in  Figure  5. 


Drag  Polar  Plot 


Figure  5:  Drag  Polar 
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The  minimum  drag  value  0.0819  occurs  at  V  =  0.9337  .  Above  this  speed,  the 
vehicle  experiences  increasing  profile  drag.  Below  this  speed,  the  aircraft  must  increase 
its  angle  of  attack  a  in  order  to  maintain  level  flight  path  angle  y.  The  required  a 
increases  more  rapidly  at  very  low  speeds,  causing  increasing  induced  drag.  Thus,  there 
are  two  values  of  velocity  for  a  given  drag  (or  thrust  level):  one  “behind  the  power 
curve”,  i.e.  at  a  lower  velocity,  and  another  one  “ahead  of  the  power  curve”.  The  effects 
of  this  phenomenon  can  be  seen  in  Chapter  7,  Figure  22,  where  V  is  reduced  to  0.9  . 
Despite  the  velocity  change,  the  required  trim  thrust  setting  ]I  remains  virtually  the  same 
because  we  move  to  an  operating  point  “behind  the  power  curve.” 

6.2.  Pitch  Dynamics  Velocity  Dependence 

We  investigate  the  V  dependence  of  the  linearized  pitch  dynamics  for  to  =  y  =  0 : 

6.2.1.  Open-loop  plant: 


a(v)= 


■2V 


K 

9^00  T74 


2 


1 

0 


b(v)= 


-2K 

-V 


Here,  the  characteristic  equation  is 


-f 


I 


+  ^  =  0 


Hence,  natural  frequency  cOn  and  damping  ratio  ^  can  be  written  as  functions  of  trim 
velocity  V : 
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4(v)=^fqC„V^+^ 


The  damping  ratio  is  minimal  at  V  =4  - - =  0.9337  and  the  minimal  damping  ratio  is 


^  =  j2KqC.  =0.0579. 


This  damping  ratio,  found  near  the  trim  velocity,  is  very  low,  and  is  the  reason  for  the 
need  of  the  damping  controller. 

Figure  6:  Open-Loop  Natural  Frequency  of  Phugoid  vs.  Velocity 

Dependence  of  Natural  Frequency  on  Velocity 
Open  Loop 
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Figure  7:  Open-Loop  Damping  Ratio  of  Phugoid  vs.  Velocity 

Dependence  of  Damping  Ratio  on  Velocity 
Open  Loop 


The  maximal  angle  of  attack  for  an  actual  F-16  is  approximately  25°  and  hence,  at 
20,000’,  we  calculate  the  stall  speed  =  0.346 .  Also,  the  approximate  value  of 
maximum  velocity  =  3 .  Plots  of  the  open-loop  natural  frequency  and  damping  ratio 

for  the  above  speed  range  are  shown  in  Figure  6  and  Figure  7. 

The  open-loop  plant  is  under-damped  for  all  values  of  velocity  of  interest.  Since 
the  region  of  interest  is  squarely  centered  in  the  region  with  lowest  damping,  the  autopilot 
must  furnish  additional  damping. 

6.2.2.  Closed-loon  Control  System 

The  proportional  controller  is  specified  by  the  gains  matrix 

^  jK^v  K„1 
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and  the  closed-loop  dynamics  are 


Act  =  A-l-BKp  - 


K,v-3 


VqCoo  +  ^  +  K-K„v 


1  +  -  2K  • 


V 


-  VK 

2  ^  ^aV 


VK 


ay 


The  characteristic  equation  is 


qC„  V  +  A  +  K .  K„,  +  i  VK„  -  i  K„, 


+  2qC„V‘K„+|5-(l  +  K„-KK„)+VK„(l  +  K„-2K.K„)-K,,.K„=0 


Recall  the  two  solutions  of  gains  option  1,  where  Kav  =  0  ,  K^y  =  0,  implying  that  the 
controller  employs  thrust  to  control  velocity,  and  AOA  to  control  flight  path  angle.  In 
this  case, 

<0.  =  ^ ^(l  -  K  •  K„ )+  qC„. V>K„  -  i ■  K„  =  V3 
^  ^{l-K  K„)+qC„  V’K„  ^ 


for  V  =  1  (trim  state).  Natural  frequency  and  damping  ratio  for  gains  option  2,  in  which 
thrust  controls  flight  path  angle  and  AOA  controls  velocity,  are  found  similarly.  Plots  of 
the  closed-loop  natural  frequency  and  damping  ratio  vs.  velocity  are  shown  in  Figures  (8) 
through  (13),  where  the  proportional  gains,  calculated  previously,  are  repeated  below: 


■K,v 

■- 0.5561 

0 

K^v 

K„y_ 

0 

1.728_ 

(Option  la) 
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omega-n 


Dependence  of  Damping  Ratio  on  Velocity 
Closed  Loop  -  GainsOption  2  -  eta: 


Figure  13:  Closed-Loop  Damping  Ratio  vs.  Velocity  (Option  2) 

Using  gain  options  la  and  lb,  the  natural  frequency  of  the  closed-loop  system  is 
fairly  constant  for  high  speeds,  but  it  changes  rapidly  in  the  low-speed  region.  The 
damping  ratio  is  greatly  increased  over  the  open-loop  case.  While  the  damping  ratio  at 

V  =  1  is  set  at  V2/2 ,  higher  speeds  cause  much  higher  damping,  and  lower  speeds  yield 
low  damping.  The  lower  damping  would  be  expected  to  cause  larger  overshoot  in  V  and 
y.  This  analysis  partially  explains  our  simulation  results  in  Chapter  7,  where  it  becomes 
apparent  that  the  operational  envelope  of  the  Mach  hold  autopilot  is  more  restricted  in  the 
low  speed  range  compared  to  high-speed  set  points. 

The  plots  for  gain  option  2  look  markedly  different.  Above  the  trim  velocity, 
natural  frequency  increases  steadily  with  velocity,  while  damping  decreases  slightly. 


59 


Regardless  of  gain  choice,  natural  frequency  and  damping  ratio  increase  rapidly  for  very 
low  velocity,  i.e.  V  <  0.5. 

Overall,  each  gain  selection  appears  to  have  improved  the  damping  ratio  across 
the  range  of  velocities  of  interest.  However,  there  are  both  areas  in  which  the  damping  is 
too  low,  and  areas  in  which  the  damping  is  higher  than  desired.  Of  particular  concern  are 
the  regions  of  low  damping,  which  could  cause  large  overshoot  by  transient  signals. 


We  investigate  the  y  dependence  of  the  linearized  pitch  dynamics  for  co  =  0 ,  V  =  1 : 


6.3.1.  Open-loot 


^n{y)=  V2  Jcos  ^  Y  +  (qC D„  +  K cos  ^  Y)sin  Y 


( 

The  damping  ratio  \  is  zero  at  y  =  a  sin 


l-^l  +  16K(qCi,^+K)^ 
4K 


=  -0.1642  =  -9.4°. 


V  J 

For  Y  <  -0.1642 ,  the  damping  ratio  ^  is  negative;  i.e.,  the  open-loop  system  is  unstable 
in  any  significant  climb.  Additionally,  since 

cos^Y  +  (qCj3^  -hKcos^ Y)sinY^0 

^  Ksin^Y  +  sin^Y“(qCD^ +K)sinY-1^0, 

does  not  exist  for  climb  angles  steeper  than  y  =  —1.351 1  =  —77.4° ;  in  other  words,  for 
Y  <  -77.4° ,  the  dynamics  are  not  oscillatory.  Plots  of  the  open-loop  natural  frequency 
and  damping  ratio  for  the  range  of —77.4°  to  90°  flight  path  angles  are  shown  in  Figure  14 
and  Figure  15. 


Dependence  of  Natural  Frequency  on  Flight  Path  Angle 
Open  Loop 


Figure  14:  Open-Loop  Natural  Frequency  vs.  Flight  Path  Angle 
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Dependence  of  Damping  Ratio  on  Flight  Path  Angle 
Open  Loop 


Figure  15:  Open -Loop  Damping  Ratio  vs.  Flight  Path  Angle 
6.3.2.  Closed-loop  Control  System 


The  proportional  controller  is  specified  by 


Kp  = 


RaV 


and  the  closed-loop  dynamics  are 


Acl  -  a  +  BKp  - 


-2(qCDo  +Kcos^Y-Kcosy-K„v)  cosy  +  K^y 

-2cosy-K„v  -sinY-K„^ 


The  characteristic  equation  is 


f  _  _  1  1  1 

+  2X  qCp^  +  K cos^  Y  +  K cos Y  •  K„v  +  - +  “Sin  Y 

+  (2qCD^  -  K^v  )(k ay  +  sin  y)  +  2(cos ^  y  +  Kpy  cos  Y  -  K  •  K cos  ^  Y  +  K  cos  ^  Y  sin  y) 
+  K„v(cosy  +  K^^  +  2K  cosy  sin  y)  =  0 


In  general, 


(0, 


'1 


{2qC^^  -  K^v  X^ay  +  sin y)+  2(cos^  y  +  cos y  -  K  ■  cos^  Y  +  K cos'  ysin y) 

+  K„v  (cos  Y  +  +  2K  cos  y  sin  y) 


11  1  - 

+Kcos^Y  +  Kcosy-K„v  +-K^ --K^v +-sin  Y 


co„ 


In  case  1,  where  Kav  =  0  ,  K^y  =  0, 


a)„  =V2^  qCo,  -^K^vj(KccY+sinY)+cos"Y-K-K^cos^y  +  Kcos^YsinY 


2_  1„  K.  1  .  . 


qCD„  +Kcos^y  +  -K^  --K^v  +-siny 


y[^qCDo-^K^vj(Kar+sinY)- 


+  cos^Y-K-K  cos^y  +  Kcos^  ysiny 


ay 


In  case  2,  where  Kccy  =  0  ,  K^v  =  0, 

a)„  =  V2^qCj3^  siny  +  cos^  y  +  cosy  +  Kcos^  ysin  y  +  K„v  (cosy  +  +  2K  cosy  sin  y) 

1  .  - 

^  qCp^ +Kcos  y  +  Kcosy-K„v +— smy 

sin  y  +  cos^  y  +  cos  y  +  K  cos^  ysin  y 

In  either  case,  setting  y-0  returns  the  expected 


0).=V3 


Plots  of  the  closed-loop  natural  frequency  and  damping  ratio  vs.  flight  path  angle  are 
shown  in  Figures  16  through  20,  where  the  proportional  gains  are 
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Kp  = 


Kp  = 


Kp  = 


1 

1 

[Kav 

K'aJ 

- 

_ 1 

[Kav 

KavJ 

- 

Kpv' 

r 

1 

<; 

7^ 

1 _ 

0.5561  0 

0  1.728 

-1.7055  0 

0  0.5787 


0  -0.9118 

31.9955  0 


(Option  la) 


(Option  lb) 


(Option  2) 

Recall  that  gains  found  by  options  (la)  and  (lb)  were  calculated  by  assuming  the 
controller  changes  thrust  to  control  velocity,  and  AOA  to  control  flight  path  angle,  while 
the  gains  found  by  option  (2)  were  calculated  by  assuming  controller  changes  AOA  to 
control  velocity,  and  thrust  to  control  flight  path  angle. 


Dependence  of  Natural  Frequency  on  Flight  Path  Angle 
Closed  Loop  -  GalnsOption  1(a) 


Figure  16;  Closed-Loop  Natural  Frequency  vs.  Flight  Path  Angle  (la) 
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Dependence  of  Damping  Ratio  on  Flight  Path  Angle 
Closed  Loop  -  GainsOption  2 


Figure  21:  Closed-Loop  Damping  Ratio  vs.  Flight  Path  Angle  (2) 
Figure  16  and  Figure  17  show  the  natural  frequency  and  damping  ratio  versus 
flight  path  angle  for  the  system  with  compensator  gains  set  as  in  option  (la).  The  peak 

Wn ,  approximately  the  desired ^3  ,  is  found  near  y  =  0 .  The  natural  frequency  drops  off 
above  and  below  this  peak,  with  a  minimum  of  0.7  in  a  vertical  climb.  Similarly,  the 
damping  ratio  is  near  its  minimum  at  y  =  0  >  increasing  above  and  below  this  point  to  a 
maximum  1.27  in  a  vertical  dive. 

Figure  18  and  Figure  19  show  the  natural  frequency  and  damping  ratio  versus 
flight  path  angle  for  the  system  with  compensator  gains  set  as  in  option  (lb).  While  the 

peak  Wn  again  is  found  in  a  slight  dive,  at  approximately  y  =  -1  the  natural  frequency 
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goes  to  zero,  while  the  damping  ratio  goes  to  infinity.  Climbs  steeper  than  this  exhibit 
unstable  dynamics. 

Figures  20  and  21  show  the  natural  frequency  and  damping  ratio  versus  flight 
path  angle  for  the  system  with  compensator  gains  set  as  in  option  (2).  The  peak  , 
approximately  the  desired ^3  ,  again  is  found  near  y  =  0 .  The  frequency  drops  off 
sharply  above  and  below  this  peak,  reaching  zero  at  approximately  -0.35  and  0.50. 
Similarly,  the  damping  ratio  is  near  its  minimum  at  y  =  0 ,  increasing  above  and  below 
this  point  toward  infinity  at  -0.35  and  0.50.  Outside  this  range,  the  system  is  unstable. 

In  summary: 

•  Gains  option  la  yields  stable  dynamics  for  all  velocities  and  flight  path  angles  of 
interest.  This  is  the  best  set  of  gains. 

•  Gains  option  lb  yields  stable  dynamics  except  for  a  small  range  of  steep  vertical 
climb  angles. 

•  Gains  option  2  yields  stable  dynamics  only  for  a  narrow  region  around  a 
horizontal  flight  path. 

6.4.  Dynamic  Reference  Signals 

We  can  deal  with  piecewise  constant  reference  signals  r,  provided  that  the  periods 
of  constancy  >  settling  time  ts  of  closed-loop  flight  control  system.  Hence,  when 
employing  a  zero-order  hold  (ZOH)  device  to  sample  the  dynamic  reference  signals  rv,  ry, 
and  To) ,  we  should  use  a  sampling  rate  for  rv,  ry,  of  at  least 

f,  =-^-  (28) 
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and,  for  ra  ,  use  a  sapling  rate  of  at  least 

fs  =—  (29) 

K 

where  t^^  is  the  settling  time  of  the  pitch  channel  variables  (the  v  and  y  states)  and  t^^  is 

the  settling  time  in  the  lateral/directional  channel  (the  heading  \ii).  Moreover, 
t  =max(t  ,t  ),  where  t,  and  t,  are  the  settling  times  for  the  V  and  y  states, 

respectively. 

Thus,  the  following  conditions  must  hold: 

.  f,=-  (30) 

precludes  the  “windup”  of  the  trim  calculation-based  control  system.  Windup  is  a 
stability  problem  with  PI  controllers  applied  to  actuators  that  can  saturate.  If  the 
proportional  channel  saturates  the  actuators,  the  integrating  channel  may  build  up  huge, 
saturating  error  output  levels.  As  the  proportional  controller  causes  the  error  to  decrease, 
the  integrating  channel  continues  to  “charge”.  The  integrator  output  does  not  begin  to 
decrease  until  the  error  changes  its  sign.  Thus,  the  actuators  can  remain  saturated  even  if 
the  output  is  at  the  desired  level. [7] 

Next,  in  order  to  achieve  good  tracking  performance,  enforce: 

•  f,  >  2fB  (31) 

This  is  Shannon’s  sampling  rate  requirement,  where  fs  is  the  bandwidth  of  the  pertinent 
reference  signal.  Thus,  condition  (31)  imposes  a  bandwidth  constraint  on  the  reference 
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signal.  Specifically,  we  require  >  2fg_ 


f  >2fg 

Sp  By 


and  f,  >  2fB_ ,  where  fg 


fg. ,  and 


f„_  are  the  bandwidths  of  the  V ,  y ,  and  (O  reference  signals,  respectively. 

■^03 

Finally,  we  require: 

•  A  small  overshoot  is  required,  so  that  the  response  of  the  dynamical  system,  e.g., 
\|/(t),  indeed  tracks  the  piecewise  constant  reference  signal  at  the  output  of  the 
ZOH  device.  This  mandates  “good”  damping  ratios,  e.g.. 


?  =—  E  =  — 

V2  ’  V2 


(32) 


Obviously,  condition  (22)  ^  output  of  ZOH  device  is  close  to  the  dynamic  reference 
signal,  e.g.,  vi/(t).  Hence,  if  the  three  conditions  (30)-(32)  are  satisfied,  good  tracking  of  a 
dynamic  reference  signal  is  achieved. 

Remark:  Note  that  both  tj  and  the  degree  of  overshoot  are  determined  by  the  linear 
control  module,  viz.,  the  gains,  Kp  and  Ki  of  the  controller.  Several  parameters  used  in 
calculating  the  gains,  i.e.,  ©n,  were  chosen  arbitrarily.  If  the  dynamic  reference  signals  to 
be  tracked  are  of  sufficiently  high  frequency  and  /  or  amplitude,  they  could  begin  to  act 
as  constraints  on  the  choices  of  these  parameters. 

6.5.  Chapter  Summary 

The  phugoid  dynamics  aircraft  model  under  consideration  is  strongly  nonlinear. 
Since  high  amplitude  maneuvers  and,  consequently,  large  excursions  in  the  state 
variables  are  contemplated,  it  is  most  instructive  to  investigate  the  linearized  dynamics’ 
dependence  on  the  state  variables.  Thus,  recall  that  the  small-signal  linear  controller 
module  is  based  on  the  linearized  plant  dynamics  model.  It  was  found  that  the  linearized 
system’s  dynamics  are  strongly  dependent  on  these  state  variables.  The  drag  polar  plot 
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showed  that  the  aircraft  flies  “behind  the  power  curve”,  i.e.,  in  a  region  where  it  must 
increase  trim  thrust  to  decrease  trim  velocity,  in  a  significant  portion  of  its  operational 
envelope.  Even  when  operated  in  closed-loop,  the  effective  values  of  cOn  and  ^  vary 
greatly  with  velocity.  The  relationship  between  these  design  parameters  and  flight  path 
angle  are  even  more  extreme;  control  options  (lb)  and  (2)  are  unstable  for  some  flight 
path  angles.  Given  the  failure  of  these  two  gain  solutions,  only  gains  option  (la)  will 
tested  with  dynamic  reference  signals  and  roll  disturbances  in  the  following  chapter. 
Finally,  an  analysis  of  requirements  for  accurate  tracking  of  dynamic  references  was 
performed. 
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7.  SIMULATION  RESULTS 


The  Phugoid  Damping  flight  control  system  is  simulated  using  the  Matlab 
Simulink  software.  Subroutines  for  the  trim  solving  equations,  the  aircraft  dynamics 
equations,  and  a  function  to  translate  the  dynamics  into  the  (x,y,z)  coordinate  system  are 
each  placed  in  one  of  Simulink’s  "S-Function"  blocks,  while  the  rest  of  the  simulation 
follows  the  structure  given  previously  in  Figure  3.  The  controller  is  first  exercised  with 
set  point  commands,  the  specifics  of  which  are  described  in  each  figure.  The  objective  is 
to  establish  the  operational  envelope  of  the  controller.  In  addition,  the  controller  is  also 
exercised  with  dynamic  reference  signals,  thus  demonstrating  phugoid  damping  control. 
For  all  simulations,  the  initial  reference  signal  (v,Y,ft))  =(1,0, 0),  while  the  set  point  or 
dynamic  commands  are  intended  to  force  the  initial  states  toward  some  desired  states.  A 
Matlab  "M-File"  prepares  the  desired  gains  option  (la,  lb,  or  2),  invokes  the  "ODE45" 
Simulink  differential  equations  solver  /  simulator,  plots  the  results,  and  adds  the 
simulation  parameters  text  descriptions  to  the  diagram. 

The  simulation  results  are  summarized  in  the  figures  below.  A  series  of 
maneuvers  types  are  simulated  with  the  intent  of  examining  the  system  performance  at 
the  limits  of  the  controller’s  operational  envelope.  High  amplitude  excursions  in 
velocity,  flight  path  angle,  and  yaw  angle  are  commanded,  followed  by  simultaneously 
commanded,  viz.,  multivariable,  set  point  changes  for  all  three  states. 

After  exercising  each  of  the  three  gain  configurations  in  this  manner,  the  roll 
disturbance  rejection  capability  is  examined  using  the  option  (la)  gains,  which  were 
found  to  yield  the  widest  operational  envelope.  Finally,  the  system's  ability  to  track  a 
dynamic  reference  input  is  determined,  thus  demonstrating  phugoid  damping  control. 
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7.1  Gain  Option  1(a) 
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Figure  22:  Decrease  Speed  to  0.9 
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Figure  24:  Decrease  Speed  to  0.41 
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Figure  25:  Decrease  Speed  to  0.40 
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Figure  30:  Vertical  Dive 
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Figure  36:  74  Degree  Heading  Change 
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Figure  37:  75  Degree  Heading  Change 
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Figure  38:  Accelerate,  Climb,  Turn 
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Figure  39:  Accelerate,  Dive,  Turn 
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Thus,  velocity  error  is  fed  back  to  modify  thrust,  while  flight  path  angle  error  is 
fed  back  to  modify  AOA. 

In  the  simulation  results,  we  see  that  high  amplitude  set  point  changes  are  well 
handled  by  the  controller. 

Figure  22  shows  the  system  behavior  when  the  reference  signal  V  is  changed  to 
0.9  .  The  new  trim  thrust  setting  p  is  the  same  as  the  old  one  (0.08),  and  the  new  trim  a 
increases  to  1.23.  As  expected,  the  aircraft  reacts  to  the  new  reference  by  decreasing 
thrust  until  the  new  trim  V  is  reached.  As  the  control  system  applies  the  new  trim  a  to 
the  aircraft  still  flying  at  the  old  V ,  the  flight  path  angle  y  increases.  The  compensator 
reacts  by  temporarily  decreasing  a  until  the  velocity  state  reaches  equilibrium.  Note  that 
the  trim  thrust  setting  p  does  not  change  because,  as  discussed  previously,  the  aircraft 
was  initially  operating  “ahead  of  the  power  curve”,  while  the  new  trim  puts  the  system 
“behind  the  power  curve.”  See,  e.g..  Figure  5. 

In  Figure  23,  the  commanded  velocity  decrease  to  0.7  is  more  significant,  and  the 
transient  is  greater  as  the  controller  dampens  the  oscillations. 

In  Figure  24,  the  new  reference  velocity  is  0.41.  In  an  effort  to  reach  this  state, 
the  aircraft  briefly  reverses  thrust,  and  then  pitches  up  nearly  vertical.  The  overshoot  in 
V  is  so  great  that  the  velocity  dips  below  V  =  0.07.  This  is  outside  the  flight  envelope  of 
the  aircraft;  the  vehicle  has  stalled.  The  equations  of  motion  are  strongly  nonlinear  near 
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V  =  0.  The  high  overshoot  is  due  to  the  low  damping  ratio  near  V  =  0.5,  as  shown 
previously  in  Figure  9.  Decreasing  the  reference  (commanded)  velocity  further,  as  shown 
in  Figure  25,  causes  the  flight  control  system  to  depart. 

Acceleration  maneuvers  are  easier  to  perform.  Figure  26  shows  the  simulation 
results  for  a  commanded  reference  V  =  2,  and  Figure  27  shows  the  reference  V  =  3 
results.  Since  damping  increases  with  velocity  in  this  region,  the  plots  show  over¬ 
damped  responses  during  acceleration  maneuvers,  as  expected.. 

The  system  response  to  commanding  a  dive  is  shown  in  Figure  28.  The  new  set 
point  Y  is  set  to  0.5,  while  V  remains  at  1.0.  The  new  trim  values  are  ]I  =  -0.4  and  a  = 
0.88.  Seeing  the  instantaneous  error  in  the  desired  flight  path  angle,  the  compensator 
“unloads”  the  aircraft,  reducing  a  to  approximately  0.  Initially,  the  aircraft  continues  to 
travel  in  a  horizontal  direction.  Given  the  reduced  (reversed)  thrust,  the  aircraft  slows  to 
below  the  reference  velocity.  In  response,  the  compensator  temporarily  adds  thrust,  until 
the  flight  path  angle  reaches  equilibrium. 

Figure  29  shows  a  steep  dive  ( y  =  1)  maneuver,  and  Figure  30  shows  a  vertical 

dive  (Y  =  ~ )  maneuver.  Despite  the  obvious  success  of  the  controller  at  achieving  the 

desired  state,  the  system  appears  to  require  large  control  usage.  The  thrust  required  in  the 
vertical  dive  is  0.95  in  reverse,  nearly  equal  to  the  weight  of  the  aircraft.  It  is 
unreasonable  to  expect  this  from  an  actual  aircraft,  even  using  speed  brakes.  Actually, 
however,  it  is  the  commanded  maneuver  that  is  extreme.  A  diving  aircraft  must  dissipate 
energy  or  accelerate.  By  maintaining  the  speed  V  =  1,  and  since  the  aerodynamic  drag  is 
insufficient,  the  speed  brakes  must  be  deployed.  In  a  dive,  increased  airspeed  would  be 
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both  expected  and  acceptable.  Given  the  increased  drag  brought  about  by  an  increase  in 
airspeed,  less  reverse  thrust  would  be  required. 

Similarly,  Figures  31-33  show  the  flight  control  system’s  response  when 
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commanded  to  perform  climbs  of  y  =  -0.5,  -1.0  and  -—,  respectively.  The  flight 

controls  system  is  able  to  perform  the  climbs  without  difficulty.  The  thrust  required  for 
the  vertical  climb  is  necessarily  greater  than  the  weight  of  the  aircraft.  This  would 
generally  require  the  use  of  afterburners. 

A  tight  right  147sec  turn  maneuver  is  shown  in  Figure  34.  Because  the  control 
variable  in  the  lateral  /  directional  channel  is  the  bank  angle  ([),  and  because  the 
longitudinal  states  are  unperturbed,  there  are  no  associated  dynamics.  This  is  not  the  case 
in  Figure  35.  In  this  simulation,  the  desired  heading  (v|io)  is  set  to  45".  The  aircraft 
initially  rolls  left,  then  levels  out  as  time  passes.  The  roll  produces  fluctuations  in  the  lift 
vector,  which  perturbs  the  longitudinal  states  V  and  y* 

Figure  36  shows  the  limit  of  the  control  system’s  ability  to  respond  to  a  heading 
change  command  without  altering  heading  rate  co.  At  a  commanded  \|/o  =  74°,  the  aircraft 
is  able  to  engage  the  commanded  heading.  However,  setting  \|/o  =  75°  causes  the  system 
to  fail,  as  shown  in  Figure  37.  This  limit  is  caused  solely  by  the  choice  of  gain  K(j, ,,, , 
chosen  in  Section  5.3.1.  Decreasing  this  gain  would  expand  the  window  of  allowable 
commanded  ij/o  ,  but  at  the  expense  of  performance  at  smaller  values  of  \|/o  .  The  other, 
more  reasonable,  approach  to  commanding  a  new  heading  is  to  command  a  heading  rate 
CO,  which  does  not  suffer  from  the  same  constraints,  being  limited  only  by  the  maximum 
turn  rate  of  the  aircraft. 
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The  final  four  figures  display  results  of  commanding  simultaneous  set  point 
changes  in  velocity,  flight  path  angle,  and  heading  rate.  The  aircraft  performs  adequately 
in  all  of  the  simulations. 
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Figure  42:  Decrease  Speed  to  0.9 
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Figure  43:  Decrease  Speed  to  0.7 
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Figure  44:  Decrease  Speed  to  0.40 
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Figure  45:  Decrease  Speed  to  0.39 
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Figure  46:  Increase  Speed  to  2 
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Figure  49:  Steep  Dive 
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Figure  50:  Vertical  Dive 
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Figure  52:  Steep  Climb 
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Figure  55:  45  Degree  Heading  Change 
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Figure  59:  Decelerate,  Dive,  Turn 


For  control  option  lb,  the  gains  are 
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Again,  velocity  error  is  fed  back  to  modify  thrust,  while  flight  path  angle  error  is 
fed  back  to  modify  AOA.  Note  that  with  these  gains,  the  thrust  setting  p  is  used  more 
heavily  to  control  perturbations,  while  the  AOA  a  is  used  less. 

The  plots  in  this  section  are  similar  to  the  plots  in  Section  7.1.  The  following 
description  serves  primarily  to  highlight  the  differences  between  the  two  control  options. 

As  seen  in  Figure  44  and  Figure  45,  the  system  is  again  unable  to  transition  to 
very  low  velocities.  The  slightly  lower  velocity  reached  by  this  set  of  gains  can  be 
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attributed  to  the  slightly  higher  damping  ratio  at  low  velocities,  when  the  control  law  (lb) 
is  used,  as  shown  previously  in  Figure  11.  This  produces  slightly  less  overshoot. 

As  shown  in  the  next  figures,  the  flight  control  system  is  able  to  perform  the 
commanded  velocity  increases  and  the  commanded  dives.  However,  Figure  52  and 
Figure  53  show  that  the  system  is  unable  to  perform  significant  climbing  maneuvers. 
Indeed,  this  was  predicted  in  Section  6.3,  where  it  was  shown  that  the  flight  control 
system  becomes  unstable  for  y  <  -1-0. 


7.3.  Gain  Option  (2) 
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Figure  60:  Decrease  Speed  to  0.925 
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Figure  65:  Dive  0.50 
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Figure  66:  Dive  0.93 
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Figure  69:  Climb  (-0.81) 

For  control  option  (2),  the  gains  are 

K,,]  r  0  -0.9118" 

|_K„v  K„J  [31.9955  0 

This  gain  set  implements  and  “opposite”  control  strategy,  as  compared  to  the  two 
previous  gain  sets.  Here,  velocity  error  is  fed  back  to  modify  AOA,  while  flight  path 
angle  error  is  fed  back  to  modify  thrust.  This  control  law  was  shown  in  Section  6.3  to 
have  a  sharply  reduced  operating  envelope.  The  simulations  verify  this.  The  control 
system  is  unable  to  decrease  speed  below  V  =  0.92  or  increase  speed  above  V  =  1.04.  It 


can  accurately  track  a  y  =  0.40  dive  reference,  but  can  not  assume  a  dive  angle  of  y  = 


0.50.  Commanding  y  =  0.93  causes  the  aircraft  to  depart  controlled  flight.  Similarly,  it 
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can  accurately  track  a  y  =  -0.30  climb  reference,  but  cannot  climb  to  y  =  -0.40,  and 
commanding  y  =  -0.81  causes  the  aircraft  to  depart. 

7.4.  Lateral  Disturbance 

In  an  effort  to  determine  the  robustness  of  the  flight  control  system,  a  lateral 
channel  disturbance  is  simulated  by  adding  a  right  roll  disturbance  to  the  aircraft  input, 
viz.,  (])  :=([)  +  (j)d  .  Two  disturbance  types  are  input:  a  constant  disturbance,  and  a  ramp. 
Each  disturbance  type  is  simulated  at  various  amplitudes.  Since  it  was  determined  in 
previous  sections  that  the  gains  calculated  by  option  (la)  are  superior,  all  simulations  in 
this  section  are  executed  using  those  gains  in  the  longitudinal  channel.  Also  as  calculated 
previously  in  section  5.4,  integral  action  is  included  in  the  lateral  channel.  No  set  point 
changes  are  commanded  during  the  disturbance  rejection  simulations. 
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Figure  70:  Roll  Disturbance  =  %I2 
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Figure  71 :  Roll  Disturbance  =  0.99-7r 
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Figure  72:  Roll  Disturbance  =  % 
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Figure  75:  Roll  Disturbance  =  2.247ft 

Figure  70  shows  the  simulation  results  for  a  high  amplitude  constant  roll 
disturbance  (l)d  =  7i/2,  applied  at  t  =  0.  This  large  value  of  (l)d  is  chosen  to  clearly  illustrate 
the  flight  control  system’s  behavior  under  a  relatively  high  disturbance. 

Initially,  the  disturbance  can  be  seen  as  a  nonzero  input  to  the  plant.  The 
aircraft’s  heading  changes  and  the  aircraft  drifts  to  the  right.  The  controller  reacts,  rolling 
the  aircraft  to  the  left  until  the  system  returns  to  the  commanded  heading;  the  aircraft’s 
path  is  displaced  in  the  lateral  direction.  Figure  71  and  Figure  72  show  the  boundary  of 
the  control  system’s  performance  capability;  if  the  roll  disturbance  inverts  the  aircraft, 
departure  may  ensue. 

Figure  73  shows  simulation  results  for  a  dynamic  roll  disturbance  (|)d  =  Ti-t.  Again, 
given  that  the  goal  of  this  chapter  is  to  determine  the  limits  of  system  performance,  the 
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simulated  disturbance  is  relatively  large.  As  the  disturbance’s  amplitude  increases,  the 
aircraft  heading  is  driven  off  the  commanded  heading.  The  compensator  eventually 
returns  the  system  to  wings-level  and  a  recovery  is  accomplished.  Note,  however,  that 
the  aircraft  settles  to  an  undesired  heading;  a  steady-state  heading  error  results.  This  is 
characteristic  of  the  ramp  response  of  a  first-order  system,  [8].  Figure  74  shows  the 
boundary  of  the  control  system’s  performance  capability,  while  Figure  75  shows  the 
departure  that  results  when  exceeding  this  boundary. 

It  is  the  integral  action  in  the  controller’s  lateral  channel  that  makes  this  recovery 
possible.  Without  it,  the  constant  disturbance  would  drive  the  system  to  an  incorrect 
heading,  while  the  ramp  disturbance  induced  error  would  be  unbounded,  causing  a  loss  of 
control  and  departure  from  controlled  flight. 

7.5.  Tracking  a  Dynamic  Reference  Signal 

In  order  to  test  the  system’s  capability  to  react  to  dynamic  reference  signals,  and 
thus  exercise  “phugoid  damping”  control,  the  previously  constant  reference  signals  were 
replaced  with  sinusoidal  ones  (in  fact,  A  cos(27tft)).  The  velocity  reference  signal 
amplitude  was  chosen  to  swing  between  the  previously  simulated  levels  of  0.41  and  3.0. 
The  flight  path  angle  reference  signal  amplitude  covered  the  full  ±n;/2  range  possible. 
The  yaw  rate  reference  signal  was  allowed  to  swing  between  ±  5,  which  was  statically 
simulated  previously,  and  is  approximately  the  aircraft’s  structural  limit.  Now,  in  the 
simulations  experiments,  the  frequency  “f  ’  of  the  sinusoidal  reference  signals  is  increased 
from  0,  which  corresponds  to  a  set  point  change  command,  to  fmax ,  where  windup 
inevitably  occurs,  viz.,  the  phugoid  damping  controller  breaks  down. 
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Figure  76:  Tracking  of  V  reference  =  1.7  +  sin(0.3  t) 
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Figure  77:  Tracking  of  V  reference  =  1.7  +  sin(0.57  t) 
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Figure  78:  Tracking  of  V  Reference  =  1.7  +  sin(0.58 1) 
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Figure  79:  Tracking  of  Gamma  Reference  =  Tsin(t) 
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Figure  80:  Tracking  of  Gamma  Reference  =  3-sin(t) 
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Figure  81:  Tracking  of  Omega  Reference  =  10-sin(t) 


Figure  76  shows  the  system  tracking  a  velocity  (V)  reference  that  changes  with  a 
frequency  f=  0.3,  which  is  nondimensionalized  in  keeping  with  previous  analysis.  Note 

that  the  dimensional  frequency  is  0.3-^  Hz  =  0.0147  Hz .  The  aircraft  tracks  the 

^0 

reference  closely.  Figure  77  shows  the  limit  of  the  system’s  ability  to  track  the  changing 
velocity  reference,  with  V’s  frequency  of  0.57.  Figure  78  shows  the  system  breaking 
down  as  it  attempts  to  track  a  reference  with  a  frequency  of  0.58.  The  higher  frequency 
of  the  reference  signal  causes  the  aircraft  to  repeatedly  reverse  direction. 

Figure  79  shows  the  system  tracking  a  flight  path  angle  (y)  reference  whose 

frequency  is  1.0  ( =  —  Hz  ==  0.05  Hz ).  The  aircraft  tracks  the  reference  closely.  Figure 
^0 

80  shows  the  limit  of  system’s  ability  to  track  the  changing  flight  path  angle  reference. 
The  Y  output  sinusoid  has  decreased  by  3  dB,  versus  a  reference  signal  y  that  changes 
with  a  frequency  of  3.0.  A  further  increase  in  this  frequency  causes  a  decrease  in  output 
amplitude,  which  signals  operation  beyond  the  bandwidth  of  the  closed-loop  system. 

Figure  81  shows  the  flight  eontrol  system  traeking  a  yaw  rate  (cOc)  varying  with  a 

frequency  of  10  ( =  10—  Hz  ~  0.5  Hz ).  The  aircraft  tracks  the  reference  closely. 

^0 

Increasing  the  frequency  of  the  yaw  rate  has  no  effect  on  the  flight  control  system.  This 
is  so  because  (i)  roll  angle  is  a  control  variable,  so  there  is  no  lag  in  the  plant’s  ability  to 
change  yaw  rate,  as  long  as  the  maximum  yaw  rate  is  attainable;  and  (ii)  maximum  yaw 
rate  is  limited  in  this  simulation  to  an  attainable  level  of  5. 
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7.5.  Chapter  Summary 


The  simulations  were  performed  using  Matlab’s  Simulink,  with  each  gain  option 
being  tested  and  the  aircraft  going  through  a  variety  of  multivariable,  high  amplitude,  and 
dynamic  maneuvers.  As  predicted,  gains  option  (la)  performed  better  than  the  others, 
demonstrating  a  wide  operational  envelope.  Deceleration  maneuvers  were  shown  to  be 
somewhat  problematic.  At  low  speed,  the  velocity  dropped  below  stall  speed.  Control 
option  (lb)  performed  adequately  in  most  maneuvers,  but  could  not  perform  steep  climbs. 
Control  option  (2)  had  a  narrow  operating  envelope  and  therefore  frequently  caused  a 
departure  from  controlled  flight.  All  results  agreed  with  the  analysis  done  previously  in 
Chapter  6. 

In  addition,  the  flight  control  system,  using  gains  calculated  using  option  (la), 
was  tested  against  lateral  disturbances.  It  withstood  a  very  large  roll  disturbance  before 
failure.  Finally,  its  ability  to  track  a  dynamic  reference  signal,  as  in  phugoid  damping 
control,  was  tested.  Even  with  the  reference  signals  having  high  amplitudes  and 
relatively  high  frequencies,  the  control  system  performed  successfully,  as  long  as  the 
(nondimensional)  reference  signal’s  frequencies  were  below  0.57  and  3.0,  for  V  and  y 
commands,  respectively. 
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8.  CONCLUSIONS 


A  novel  phugoid  damping  control  design  methodology  is  developed,  based  on  the 
use  of  wind  axes  and  a  point-mass  aircraft  model.  A  multivariable  set  point  controller  is 
designed  which  consists  of:  (i)  A  trim  calculation-based  nonlinear  feed-forward  control 
computer;  thus,  given  a  commanded  new  trim  state  (air  speed,  flight  path  angle,  and  yaw 
rate),  the  required  trim  thrust  setting  and  trim  angle  of  attack,  roll  angle,  and  sideslip 
angle  inputs  are  solved  for,  and,  (ii)  a  small  signal  linear  feedback  regulator;  the 
equations  of  motion  linearized  about  the  trim  condition  of  wings  level,  constant  altitude 
flight,  which  simplifies  the  dynamics  to  allow  separation  between  the  lateral  and 
longitudinal  control  channels,  are  used,  and  a  small-signal  linear  multivariable  regulator 
is  designed.  The  linear  compensator  also  entails  integral  action.  Thus,  the  controller 
consists  of  a  strongly  nonlinear  feed-forward  module  and  a  linear  small  signal 
compensator.  The  novel  proposed  multivariable  set  point  controller  encompasses  full 
three-axes  autopilot  functions.  The  command  signals  are  airspeed,  flight  path  angle,  and 
heading  angle  or  heading  rate.  Moreover,  this  controller  is  used  as  a  tracking  controller, 
a.k.a.  a  “phugoid  damping”  controller,  provided  that  the  bandwidth  of  the  command 
signal  is  substantially  less  than  the  bandwidth  of  the  closed  loop  flight  control  system. 
Robust  and  flexible,  the  flight  control  system  exhibits  the  capability  to  independently 
control  either  bank  angle  or  sideslip  angle,  allowing  the  air  vehicle  to  perform  weapons- 
pointing  tasks.  The  phugoid  damping  controller’s  performance  is  examined  in  extensive 
simulations  and  its  wide  operational  envelope  is  demonstrated.  It  is  shown  that  the 
controller  can  accept  high  amplitude  commands.  Thus,  the  speed  can  be  commanded  to 
change  between  V  =  0.41  and  V  =  3.0 ,  the  flight  path  angle  can  be  varied  anywhere  in 
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_  71 

the  full  range  of  y  =  ±— .  and  heading  changes  can  be  offset  ±  74° ,  while  heading  rate 
commands  are  limited  only  by  the  maximum  turn  rate  of  the  aircraft. 
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